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MODULUS  DATA  INPUT  FOR  THE  COMPUTER 


The  relaxation  modulus  was  represented  by  a  sixteen-term  Prony 
series.  This  relation  for  the  tensile  modulus,  E(t),  is 


m  -3  t 

E(t)  =  A  +  £  A,  e 
0  i=l  1 


(A-l) 


where  Aq,  A^,  and  3^  are  constants, 

For  an  -ncompressible  material  the  relaxation  modulus  in  shear,  p 
is  given  by 


h(t) 


E(t) 

3 


m 

+  £ 
i=l 


(A-2) 


The  viscoelastic  stress  analyses  requires  the  shear  modulus 
representation. 

The  Prony  series  constants  employed  in  the  present  parameter 
study  are  given  in  Table  A-l,  for  the  CTPB  propellant,  and  in  Table  A-2, 
for  the  HTPB  propellant. 

To  complete  the  viscoelastic  description  the  time-temperature 
shift  function,  a  ,  is  required.  Hence,  also  listed  in  Tables  A-l 
and  A-2  are  the  logarithims  of  a  at  12  different  temperatures  from 
-100c  to  +200°F,  for  the  two  referenced  propellants. 
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MODULUS  INPUT  FOR  THE  HTPB  PROPFI.LANT 


Prony  Series  Parameters 


Log  (t/a  ) 
(min) 

E 

psi 

(hr 

■h 

Ai 

A./3 

i 

q 

80 

26.67 

0 

-8 

10000 

3  x 

108 

107 

4376.65 

1458.9 

1 

-7 

7000 

3  x 

3697.78 

1232.6 

2 

-6 

4000 

3  x 

106 

2056.49 

685.50 

3 

-5 

2400 

3  x 

105 

896,865 

298.96 

4 

-4 

1600 

3  x 

104 

665.701 

221.90 

5 

-3 

1050 

3  x 

103 

384.337 

128.11 

6 

-2 

720 

3  x 

10, 

250.028 

83.343 

7 

-1 

500 

3  x 

101 

180.589 

60.196 

8 

0 

350 

3  x 

100 

10-1 

10-2 

106.550 

35.517 

9 

1 

2 

260 

205 

3  x 
3  x 

67.0877 

35.8845 

22.363 

11.962 

10 

21 

3 

170 

3  x 

10-3 

39.3243 

11.441 

12 

4 

140 

3  x 

24.4644 

8.155 

13 

5 

120 

3  x 

10  J 

15.1300 

5.043 

14 

6 

110 

3  x 

10  J 

-2,4636 

-0.8212 

15 

1 

100 

3  x 

10  6 

33.0018 

11.001 

16 

Time-Temperature  Shift  Function 


Temp, ,  °F 


L°810aT 


-75 

6.90 

-50 

4.93 

-25 

3.41 

0 

2,35 

25 

1.48 

50 

0.72 

77 

0 

100 

-0,46 

125 

-0.90 

150 

-1.24 

170 

-1.47 

A- 3 
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PARAMETER  STUDY  FOR  HISTORY  1 


The  results  of  a  series  of  one-dimensional,  thermoviscoelai-ic 
analyses  are  presented  graphically  in  this  appendix.  All  of  the  analyses 
were  based  upon  the  environmental  temperature  history  described  In  the  text 
and  shown  separately  in  each  figure.  The  data  for  the  CTPB  propellant  are 
presented  in  Figures  B-l  to  B-6,  while  those  for  the  HTPB  propellant  are 
given  in  Figures  B-7  to  B-12. 

No  analyses  of  the  results  are  made  here. 

A.  CTPB  PROPELLANT 

The  following  giaphs  give  the  thermoviscoelastic  solutions  for 
these  grain  designs. 


B-l. 


Strain 
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PARAMETER  STUDY  OF  INNER- BORE  HOOP  STRAINS 


CTFB  PROPELLANT  -  HISTORY  1 
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Figure  B-l 
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PARAMETER  STUDY  OF  INNER-BORE  HOOP  STRAINS 
CTFB  PROPELLANT  -  HISTORY  1 
(B  =  8  in.) 
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Figure  B-2 
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PARAMETER  STUDY  OF  INNER-BORE  HOOP  STRESS  SOLUTIONS 
CTPB  PROPELLANT  -  HISTORY  1 
(B  =  4  in.) 
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Figure  B-3 
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PARAMETER  STUDY  OF  INNER-BORE  HOOP  STRESS  SOLUTIONS 
CTPB  PROPELLANT  -  HISTORY  1 
(B  =  &  in.) 
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Figure  B-4 
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B.  HTPB  PROPELLANT 

The  following  graphs  give  the  thermoviscoelastic  solutions 
for  these  grain  designs. 


B-8 
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HTPB  PROPELLANT  -  HISTORY  1 
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Figure  B-7 
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PARAMETER  STUDY  FOR  HISTORY  2 


The  results  of  a  series  of  one-dimensional,  therraoviscoelastic 
analyses  are  presented  graphically  in  this  appendix.  All  of  the  results 
were  based  upon  the  environmental  temperature  history  described  in  the  e 
text  and  shown  separately  in  each  figure.  The  data  for  the  CTPB  and  1ITPB 
propellants  are  presented  separately. 

No  analyses  of  the  results  are  made  here. 

A.  CTPB  PROPELLANT 

The  following  graphs  give  the  thermoviscoelastic  solutions 
for  these  grain  designs. 
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PARAMETER  STUDY  OF  RADIAL  BOND  STRESS  SOLUTIONS 
CTPB  PROPELLANT  - ^HISTORY  2 

(B  =  4  in.) 
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B.  HTPB  PROPELLANT 

The  following  graphs  give  the  thermoviscoelastic  solutions 
for  these  grain  designs. 
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PARAMETER  STUDY  OF  INNER-BORE  HOOP  STRESS  SOLUTIONS 
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(B  =  4  in.) 


Aj  B/A 


2.0 


Stress 


Aerojet  Solid  Propulsion  Company 
Report  1341-26F 


Appendix  C 


PARAMETER  STUDY  OF  INNER-BORE  HOOP  STRESS  SOLUTIONS 
Hl’PB  PROPELLANT  -  HISTORY  2 
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Figure  C-10 
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A  NEW  NORMALIZED  RELATION  FOR  THE  RELAXATION  MODULUS 


A.  NORMALIZATION  OF  THE  SHEAR  MODULUS 

Conventional  viscoelastic  stress  analyses  involve  the  relaxation 
modulus  in  shear,  u(t).  The  basic  relations  for  this  modulus  require 
difficult  to  perform  experiments,  the  results  of  which  are  not  always 
satisfactory.  It  was  desirable,  therefore,  to  devise  a  simpler  to  use 
expression  for  this  modulus.  This  was  done  starting  with  the  relation 
now  used , 


u(t)  =  p  +  -  u  )  /  he  t//T  d£m  (D-l) 

e  8  e 

where  p 

g 

Ue 
h 

t 
T 


is  the  glassy  shear  modulus 
is  the  equilibrium  shear  modulus 

represents  a  continuous  distribution  of  relaxation  times 
(normal! zed) 

is  the  time  of  the  test 

is  a  relaxation  time 


Since  u  is  difficult  to  evaluate  experimentally  an  attractive 
substitute  wa§  sought  and  found. 

We  solved  for  p(t)  at  some  specific  time,  like  one  minute,  to 
obtain  p(l). 


u(l)  «  Pe  +  <pg 


V 


/  he 


-l/i 


d&m 


(D-2) 


Letting 

C  -  f  he  d?nr 

—  c*> 


(D-3) 


where  C  is  a  constant,  we  solve  for  p  -  p  in  Equation  (D-2)  using  Equation 
(D-3)  and  inserting  the  result  into  Equation  (D-l)  gives 

-l/idim 

’.:  ( t )  =  u  +  (p(l)  -  u  )  ~  (D-4) 

e  e  C 

— oo 


D-I 
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For  practical  experimental  purposes  the  constant  C  can  be  combined 
with  h  to  give  h',  a  quantity  which  is  experimentally  identified  in  the 
same  manner  as  h. 

Thus,  the  new  relation  becomes 

p(t)  =  vig  +  (p(l)  -  pg)  /  h'e  d£m  (D-5) 


In  engineering  practice  is  is  unnecessary  to  evaluate  h'.  Instead, 
a  graphical  plot  of  p(t)  versus  time  is  usually  sufficient.  When  normalized 
results  are  required  a  plot  of  p(t)  -  ue  versus  the  time  is  equivalent 

ToT^lT 

00  —  t/T 

to  a  plot  of  /  h'e  d£nt  versus  the  time. 

r  CO 

Obviously,  a  broad  range  of  relaxation  curves  can  be  obtained  for 
a  given  distribution  of  relaxation  times,  h1. 

B,  NORMALIZATION  OF  THE  PRONY  SERIES 

The  new  relation  permits  a  normalization  of  the  Prony  Series  as 
well.  This  relation  in  its  usual  form  is  given  as 


m  -6  t 

p(t)  =  a  +  T.  a  e  m  (1)~6) 

o  -  m 


where  the  a  and  B  are  constants 
m  m 

«q  is  the  equilibrium  relaxation  modulus 


We  can  normalize  Equation  (D-6)  to  give  a  form  similar  to  that  of 
Equation  (D— 5 ) .  First  we  replace  a  by  the  equivalent  term  \i  then  normalize 
the  constants  am,  as  shown  below,  to  give 

m  -6  t 

p(t)  =  Pe  +  (y(l)  -  Pg)  £  a'ni  e  m  (D-7) 

m=l 


where 


a  /  (p(1)  -  u  ) 
m  e 


(0-8) 


D-2 
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Equation  (D-7)  forms  the  engineering  basis  of  our  normalization  procedures. 
This  normalization  method  defines  the  uiulus  in  terms  of  two  easily  determined 
parameters,  p  and  n(l)  -  p  .  These  same  parameters  can  be  used  to  normalize 
the  stress  ant  strain  data  from  our  engineering  analyses. 
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INCREMENTAL  ANALYSIS  PROCEDURE 


In  an  attempt  Lo  minimize  accumulated  numerical  errors  in  our 
linear  viscoelastic  analyses,  Dr,  Herrmann  developed  this  new  approach. 

The  "total  stress  analysis"  was  replaced  by  an  "incremental  analysis 
procedure";  i.e.,  instead  of  solving  fcr  che  total  stress  and  strain  in 
the  propellant  for  a  given  point  in  time,  one  could  solve  instead  for  the 
incremental  changes  in  stress  and  strain.  The  incremental  equations  are  very 
similar  to  those  previously  reported  for  the  "total  analysis"  with  the 
following  exceptions  (the  equations  referred  to  by  number  are  reported  in 
Reference  (E-l). 

Consider  first,  Equ ttion  (13), 


2pN  eijN  +  LijN 


(13) 


The  quantities  S. .  and  e, .  in  Equation  (13)  need  to  be  interpreted 
1JN  1JN 

as  the  incremental  changes  in  stress  and  strain  during  time  step  N  (i.e., 

S..  =  2u..Ae..  +  L.  .  where  S..  =  S.,  +  AS..  ,  etc.).  Equation 

l-’N  N  1JN  1JN  lJN-l  1JN 


LIJX  ~  2  ^XilN  ~  ^PN  "  ao)  eilN-lJ 


is  replaced  by 


E-l 
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Equation  (19) 


M 

XijN  =J'_X  %  CijNm 


(19) 


is  eliminated  and  Equation  (24) 


-6  AF,,, 

m  N 


CijNm  S 


ICijN-l,m  +  (eijN~l  "  eijN-2*  Vl.m1 


(24) 


becomes 


-(3  A?..  , 

C.  .  -  e  “  “-1  C 

i;lH« 


.  ,  +  Ae .  J.,  , 

^N-l.m  ljN-l  N“1,m 


(24) 


Finally,  Equations  (1),  (2)  and  (3) 


Tij  Sij  *  6ij  a>  0  =  33  Tii 


(1) 


eij  eij  +  6ij  3  ’  0  "  eii 


(2) 


a  =  K  (9  -  3a AT) 


(3) 


where  K  =  Bulk  modulus 

a  =  Coefficient  of  linear  thermal  expansion 
AT  =  T  (x,  t)  -  T 

T&  =  Initial  stress  free  reference  temperature 
are  replaced  by  their  corresponding  incremental  forms  (Note:  AT^.  =  T  -  ^ ) . 


E-2 
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PROXY  SERIES  CURVE  FIT  ANALYSIS 


The  shear  relaxation  modulus  for  most  solid  propellants  has  been  found 
to  fit  the  following  series  function  with  relatively  few  terms: 

n 


0(t)  =  A  +  Z  A.e 
°  1-1  1 


(F-l) 


Equation  (F-l)is  a  "Prony"  series  with  two  unknown  coefficients  A.  and  (i. 
The  method  of  collocation  is  used  to  find  these  coefficients:  1 


Let  6 


1 

*1 


( F— 2 ) 


Substituion  of  Equation  (F— 2)  into  Equation  (F-l)  gives 


_ t 

2t\ 


0  =  A  +1  A  e 
°  i-1  1 


(F-3) 


Now,  choose  n  points  for  the  evaluation  of  0 

t. 

n  _ JL 

0.  =  A  +  1  A.e  2t, 

1  0.-1  l 

J  1=1 


(J  -  1  >n) 


(F-A) 


Equations  (F-4)  are  sufficient  to  solve  for  the  A^.  In  matrix  notation! 


[E'j  {A.}  =  {0^  -  Ao) 


(F— 5) 


where , 


E. . 
Ji 


e 


2ti 


Equation  (F-5)  is  readily  solved  for  A^; 

fA.l  =  lEf]{0.  -  Ao) 


( l’“6 ) 


F-l 
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Equation  (F-6)  has  been  programmed  for  computer  solution  from  a  times 
burning  terminal  in  the  BASIC  language,  A  listing  of  the  program  is  given 
on  Page  F-3.  The  order  of  data  input  is  given  below: 


1.  n 

2.  A 

o 

3.  tj  (j  =  1,  n) 

4.  0j  (j  -  1,  n) 

Data  statements  600  to  9000  may  be  used  for  data.  Sample  data  statements 
are  shown  below: 

600  Data  7,100 

610  Data  IE-4,  IE-3,  IE-2,  IE-1,  1,  10,  100 
620  Data  3000,  1800,  1000,  620,  410,  320,  280 
A  sample  run  with  this  data  is  shown  on  Page  F-4. 
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PROGRAM  LISTING 


1  00 
I  JO 
1  2  0 
1  3C 
1  40 
1  SO 
1  60 
1  70 
1  8G 

1  90 

2  00 
0  to 
2  20 
2  3C 
2  40 
2  SO 
0  60 
2  70 
2  3C 

2  90 

3  00 
:<  10 
0  21' 
3  30 

40 
J  50 

3  6( 
:■  7o 
1  SC 

>0 

4  00 
4  I  Cl 


PFM  "PPTNY  SFPIFS  CURVF  FIT'* 

blM  T  (  20,  I  > ,  M  20,  1  > ,  AC  20,  20)  ,  CC  20,  2  0) 

Mv  Ml,  20) ,  C  (  1,20) 

PEAL  N  ,  fe  1 


PEAL  lCN,i) 

P  =  ZF.PC  1  ,N  > 

C  =  ZF.PC  i  »N  > 

F  =  TFM(T) 

INPUT  TIMES  AND  MODULI 
PRINT  F  5 
A  =  ZEPCN,N) 

J  =  2  TO  N 

b  =  -TC  J,  1  >/T<  t,  1  ) 


ACJ,J)  = 
0 


MAT 
MAT 
MAT  C 
MAI  F 
FPlN  1 
MAT 
MAT 
F  :1  F 
LET 
LIT 

IF  b  <  -20 
LET  A( 1 . J) 

LET  AC J, 1 ) 

O  T3  300 
LF  T  AC  1  , J) 

LF  T  AC J,  1  > 
NOT  J 
LET  AC  1  .  1  ) 

LET  K  l  =  M 
RIP  1=2  13  K 1 
LF  r  L  =  I 
LFT  K2  =  K 1  - 

FTP  J  =  ir  TM  K2 
LFT  L  =  L  «•  1 


0.  606S3 
THEN  280 
EXPCO.  S/P) 
F  X  PC  0.  S*E  ) 


0 

0 

.  606S3 


LET  A(  j,|.)  =  A<  J,  l  ) 
LFT  ACL, J)  =  AC  1 ,  1) 
Nr  XT  J 
'J  t  K  T  I 


4  20  MAT  C  =  ZF.RCN.N) 

<.30  MAT  C  =  IN  VC  A) 

4  40  MAT  PEAL  P.CN,  1  ) 

'  SO  MAT  (  =  1RNCE  ) 

4  60  MAT  PP  IN  T  C  ) 

4  70  F3P  .1=1  TM  N 
4  S  o  i.FT  E  C  J,  1)  =  P  C  J,  1  )  -  FI 
4  9  0  N  F  X  T  J 


SCO  MAI  T  =  OF 
S  10  PP1NT  "  MLUT  I MN  " 
S  20  vat  F  =  TRN(T) 

S  JO  vf.i  p P  1  N  T  F  J 


4  4('  M  7  1 
-i  ■  •  >  -S’. 
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SAMPLE  RUN 


P  R8N  9:57  SF  FBI  02/13/70 

I  iN?UT  TIMES  AN  1/  MODULI 

.  0001  .  001  .01  .1  110  100 

3000  18  00  1000  620  410  320  280 

S0LUTI3N 

1370.27  1031.23  464.  66  247.  646  184.314  -106.  615  297.955 

BUT  3F  DATA  IN  440 
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INCLUSION  OF  NON-ZERO  THICKNESS  STRESSES 
IN  PLANE  STRESS  ANALYSES 


Th'i  generalized  plane  stress  option  was  modified  so  that  the  stress 
throughout  the  thickness  of  the  body,  may  be  specified  as  a  constant  rather 
than  zero.  The  values  of  the  thickness  stress  constitutes  one  of  the 
input  parameters.  The  inclusion  of  a  non-zero  thickness  stress  required 
the  modification  of  the  governing  variational  equation. 

It  was  determined  that  the  appropriate  form  of  the  variational 
function,  for  a  nearly  incompressible  material*’  '  is  (for  T  =  T^) : 


fn  ■  m  "T  [<4Ex>2  +  %  /  -  +  *  M» 

N  N  «  N  N 


+  (u„  AIL  -  A0../3)  (Ae  +  Ae  )  +  L  Ac  .  +  L  Ae 
*  ^  N  *8  XXM  XN  + 

.2 


+  L  Ay 
xyN  XyN 


2S 


-  2  -  Af5N 

(AHj,)^  +  AHjj  (rf  -  3pN  a  AT^) 


-  AF  Au  AF  Au  }  dxdy  -  /AuM  • 


0-1 


A (applied  boundary  load)  ds 
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where 


2"n 

s=  - 

xx,,  3 


AB. 


(2 Ac  -  At  )  +  p  AH  +  L 


'N  'N 


At 


AAli. 


‘  (At  +  At  ) 
3  XN  yN 


N  N  xx 

h 

+  t2~  (Ao  -  L  )  + 

3un  N  Z,n 


3  “  ATN 


All..  -  2u  ,  a  AT  +  (Ao  -  L  ) 
N  N  N  N  zz., 

N 


S  = 


‘V1  " 


4u, 


3(1  + 


3K 


The  symbol  Aa^  denotes  the  incremental  change  in  the  specified  thickness 

stress,  L  denotes  the  history  term  associated  with  c  ,  the  other 
ZZN  Z 

symbols  have  their  usual  meanings. 
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A  COMPUTER  PROGRAM  FOR  VISCOELASTIC  SOLIDS  OF  REVOLUTION 
SUBJECTED  TO  TIME-VARYING  THERMAL  AND  MECHANICAL  LOAD  ENVIRONMENTS 

-  VERSION  2.1  - 


I .  INTRODUCTION 


The  purpose  of  the  computer  code  described  in  this  section  is  to  perform 
viscoelastic  stress  analyses.  The  analyses  is  applicable  to  arbitrary  revolved 
solids  and  plane  structures  subjected  to  loads  of  mechanical  or  thermal  origin. 
The  program  is  segmented  into  two  (2)  phases:  (1)  Transient  Heat  Transfer 
Analysis;  and  (2)  Viscoelastic  Stress  Analysis.  The  purpose  of  the  heat 
transfer  phase  is  to  generate  temperature  distributions  in  the  body  as  a 
function  of  time  which  are  used  subsequently  in  performing  the  stress 
analysis.  The  two  phases  of  the  program  can  be  used  in  sequence  within  a 
given  job  or  each  phase  can  be  used  separately. 

The  method  of  solution  employs  the  finite  element  procedure  for  solving 
the  spatial  problems  (heat  conduction  and  stress  analysis)  and  time  marching 
techniques  to  evaluate  temperatures/stresses  at  successive  points  in  time. 

The  transient  heat  transfer  problem  is  solved  using  the  procedure  developed  by 
Wilson  and  Niekell*.  Knowing  thermal  and  mechanical  loads  as  a  function  of 
time  and  having  available  the  viscoelastic  properties  of  the  material (s), 
a  set  of  equivalent  elastic  parameters  is  defined  for  a  particular  point  in 
time;  the  equivalent  elastic  problem  is  posed  using  the  procedure  given  by 
Herrmann  and  Pet* "son**.  An  elastic  stress  analysis  must  be  performed  at 
each  point  in  tine  for  which  the  viscoelastic  response  of  the  body  is 
required.  The  stress  analysis  problem  is  solved  using  a  finite  element 
method  given  by  Herrmann.*** 

II.  APPLICATIONS 


The  principal  application  of  the  program  is  to  stress  analysis  of 
solid  propellant  grains  maintained  in  time  varying  temperature  environments. 

A  typical  application  might  be  the  grain  stress  analysis  of  a  motor  system 
subjected  to  thermal  cycling.  Mechanical  loads  may  also  be  applied  to  the  body 
either  isothermally  or  with  simultaneous  time  varying  temperature  changes.  The 
kinds  of  mechanical  loading  considered  in  the  analysis  include  surface  pressures, 
body  forces  (spin  and/or  axial  accelerations),  concentrated  nodal  forces,  and 
specified  nodal  displacements;  all  mechanical  loads  can  vary  with  time  in 
accordance  with  .'ser-supplied  table. 


*  Wilson,  E.  L. ,  and  Niekell,  R.  E. ,  "Application  of  the  Finite  Element 
Method  to  Heat  Conduction  Analysis",  Nuclear  Engineering  and  Design  4 
(1966),  p.  276-286.  North  Holland  Publishing  Company,  Amsterdam. 

**  Herrmann,  L.  R. ,  and  Peterson,  F.  E. ,  "A  Numerical  Procedure  for  Visco¬ 
elastic  Stress  Analysis",  CPIA  -  7th  Meeting,  Working  Croup  on  Mechanical 
Behavior  (ICRPG),  November  1968. 

***  Herrmann,  L.  R. ,  "Elasticity  Equations  for  Incompressible  Materials  by 
A  Variational  Theorem",  Journal  of  the  AIAA,  p.  1896-1900,  October  1965. 
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A.  PHYSICAL  PROPERTIES 

The  program  can  be  used  to  perform  simple  elastic  solutions, 
time-dependent  elastic  or  thermoelastic  analyses,  or  thermoviscoelastie  analyses. 
A  viscoelastic  analysis  requires  a  complete  material  iroperty  characterization 
including  the  master  relaxation  curve  and  shift  function  for  each  time 
dependent  material,  bulk  modulus,  expansion  coeffici:nt,  density,  etc.  Tabic 
H-  1  illustrates  the  thermal  and  mechanical  propert  es  data  required  to  perform 
the  thermoviscoelastie  stress  analysis  of  a  bipropeJ lant  grain  with  a  separate 
viscoelastic  liner  and  an  elastic  case. 

Shift  function  data  is  accepted  by  the  program  in  the  form  of  a 
table  of  log  aT  versus  T,  (°F);  shift  factors  at  temperatures  other  than 
those  supplied  in  the  input  table  are  determined  using  linear  interpolation 
between  given  points. 

Relaxation  data  must  be  input  as  the  coefficients  (A^,  3^) in  a 
Prony  Series  fit  to  the  experimental  data.  The  relaxation  behavior  in  shear 
must  be  expressed  as  an  exponential  series. 

n  e“6it 

<|>(t)  -  A  4-  £  A,  c  (H-l) 

°  i-i  1 


where  A  is  the  shear  equilibrium  modulus  and  (Aj.,  (^1  are  founu  from  curve 
fitting^  calculations  based  on  experimental  data.  Each  viscoelastic  material 
muse  have  Its  relaxation  behavior  expressed  in  a  separate  series  expansion. 
Elastic  materials  have  no  terms  in  the  series  other  than  A  .  If  there  are  no 
terms  in  the  series  expansion  (M  =  0)  the  program  will  not  read  a  shift 
function  table. 

The  special  problem  of  bulk  rapid  pressurization  of  a  propellant 
grain  where  a  pressure  shift  function  a  is  required  is  handled  in  a 
different  manner.  This  type  of  problemPcan  only  be  run  isothormally  with  nc 
superimposed  thermal  loads.  In  this  case  the  shift  function  input  values  are 
interpreted  as  the  product  a  a  as  a  function  pressure  where  a  is  a  constant 
for  the  temperature  under  consideration.  ^ 


*  On*,  recommended  curve  fitting  procedure  is  a  "collocation  method"  originated 
by  Schapery  and  summarized  in  the  1CRPG  Solid  Propellant  Mechanical  Behavior 
Manual  starting  at  Section  2.2-4  (June  1963).  The  method  involves  assuming 
values  for  the  constants  0^  and  solving  a  set  of  linear  simultaneous  equations 
for  the  constants  A^.  A  time-share  routine  called  "PRONY'"  has  been  programmed 
to  perform  the  calculations  at  Aerojet. 
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PHYSICAL  PROPERTY  INFORMATION  REQUIRED 
FOR  A  THERMOVISCOELASTIC  STRESS  ANALYSIS 


PROBLEM  IDENTIFICATION: 


Propellant  No.  1  - 
Propellant  No.  2  - 
Liner 
Case 


Physical  Property 

Propellant 
No.  1 

Propellant 

No.  2 

Liner 

Case 

Coefficient  of  Linear 
Thermal  Expansion  (°F-  ) 

Density  (lb-in.  ") 

Specific  Heat  . 

Capacity  (btu-lb"  °F~  ) 

Thermal  Conductivity, 
(Btu-in.“^  hr-^  °F--) 

Bulk  Modulus  ( lb-in. ~  ) 

(Table  or  curve  as  a  function  of 
time) 

,  -  .  .  _  - 

Shift  Factor  (a^  vs.  °F) 

(Table  or  curve  as  a  function  of 
temperature) 

Elastic  Shear 

Modulus  ( lb-in. 

— 
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III.  DATA  INPUT  DESCRIPTION 

This  section  supplies  information  necessary  for  the  preparation  of  data 
input  cards.  The  input  sequence  is  separated  into  four  major  groups: 


1.  Grid  Definition 

2.  Solution  Time  and  Temperature  Information 

3.  Transient  Heat  Transfer  Solution  Data 

4.  Stress  Analysis  Information 


A.  NOMENCLATURE 

The  abbreviation  "cc"  used  below  stands  for  "card  columns".  The. 
variable  names  assigned  to  the  various  parameters  used  by  the  program  are 
given  below  in  upper  case  letters;  for  example,  "NNP"  stands  for  the  total 
number  of  nodes  in  the  finite  element  mesh.  All  variables  starting  with 
any  of  the  letters  I,  J,  K,  L,  M,  N  are  to  be  input  to  the  program  as 
integers  (i.e.,  without  a  decimal  point).  All  integers  are  to  be  packed 
to  the  right  of  the  field  specified  by  the  "cc"  numbers.  Any  variable  whose 
first  letter  is  not  an  I,  J,  N  is  a  real  number  requiring  a  decimal; 

"R(N) for  example,  is  the  radius  of  the  "N-th"  nodal  point  (entered  in  cc 
6-15).  If  R(N)  =  13.45,  then  the  number  "13.45"  can  be  placed  anywhere  in  the 
field:  cc  6-15.  Real  numbers  can  also  be  input  in  "E"  formal:;  13.45  could  be 
entered  as  1.345E1,  .1345E2,  1345. E-2,  etc.,  providing  the  set  of  characters 
is  packed  to  the  right  of  the  cc  field. 

Whenever  applicable,  units  of  the  variables  are  stated  in 
symbolic  notations.  The  symbols  used  below  are  defined  as  follows: 
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(?)  =■  Force  units 

(L)  =  Length  units 

(T)  --  Time  units 

(°F)  =  Temperature  units 

(R)  =  Radian 

(Btu)  =  Thermal  Heat  Flow  Units 

Thus,  the  quantity  (F)  (L) "  written  after  the  shear  modulus  means  *psi!*  if  pounds  and 
inches  are  the  units  chosen  by  the  user.  There  are  no  units,  conversion  factors, 
etc.,  built  into  the  program;  thus,  once  a  set  of  units  is  chosen  it  must  be  used 
consistently  throughout  the  analysis. 

• 

B.  SEQUENCE  OF  OPERATIONS 

The  flow  of  program  execution  is  controlled  by  three  (3)  user-supplied 

variables : 


IPF  =*  Plot  Control  Flag 

NT3M  =  Temperature  Information  Flag 

J$B  =  Job  Control  Flag 


Table  H-2  contains  values  which  can  be  assigned  to  these  variables 
showing  what  operation  results  from  a  particular  specification.  Certain  combinations 
are  not  possible.  For  example,  if  XPF<0,  the  program  reads  mesh  data,  prepares  a 
plot  and  returns  control  to  that  portion  of  the  program  which  looks  for  the  next 
job;  thus,  the  variables  "NTEM"  and  •JOB*  cannot  be  specified.  *JOB"  has  no  meaning 
unless  IPF  >  0  and  NTEM  =  0. 

The  option  NTEM  =  1  is  impractical  for  real  problems  because  of  the  amount, 
of  card  data  involved;  this  option  is  useful  when  solving  * check  cases*  with  temperature 
distributions  generated  from  an  analytical  expression  or  formula.  The  NTEM  =  2  option  saves 
re-runing  the  temperature  problem  if  only  the  mechanical  properties  or  mechanical  loads 
change.  NTEM  =  3  is  used  for  the  isothermal  problem  in  which  the  loading  is  mechanical 
in  origin.  Viscoelastic  materials  exhibit  temperature  dependence,  so  material  temperatures 
must  be  defined  even  though  there  are  no  driving  thermal  strains  in  the  body  to  be 
analyzed. 


The  va  riable  "J$B"  (which  is  only  defined  if  NTEM  =  0)  controls  what  the 
program  does  with  the  results  of  the  heat  transfer  solution  and  where  the  program 
goes  after  the  heat  transfer  calculations,  if  J$B  =  0,  element  temperatures  are  saved 
temporarily  for  use  in  the  stress  analysis  to  follow.  J0B  =  1  results  in  the  same 
operation  as  J$B  =  0,  but  in  addition  the  element  temperatures  are  saved  for  use  in  a 
later  analysis  (or  series  of  analyses).  If  J$B  =  2,  the  temperatures  are  saved.;(and 
printed  at  user  specified  time  intervals)  on  a  tape  for  use  at  some  other  time;  at  this 
point  the  program  is  through  with  this  job  and  looks  for  another.  J0B  =  3  allows  the 
user  to  run  the  program  solely  as  a  heat  transfer  analysis1;. 
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EXECUTION  CONTROL  VARIABLES 


FnAG 

VALUE 

OPERATION  PERFORMED 

li’F 

<  0 

Read  grid  data  only,  plot.  .rid  and  :.ton. 

=  0 

Read  all  data,  execute  the  jt without  t  plot. 

>  0 

Read  all  data,  execute  the  job  with  a  plot. 

NTEM 

=  0 

Calculate  the  I'lemant  temperature:;  using  the 
heat  transfer  analysis. 

=  1 

Read  the  element  temperatures  from  card  input. 

a  2 

Read  the  element  tempera tures  from  a  tape  rivaled 
during  a  previous  run. 

«  3 

.  . 

Read  the  element  temperatures  from  card  inpui. 
and  use  the  sumo  distribution  for  all  solution 
time  points  (isothermal  response). 

=  0 

Run  the  heat  transfer  problem  and  use  the  iv/nit. 
to  perform  the  stress  analysis. 

=  1 

Run  the  heat  tran  for  problem,  save  Uio  •  i«-m.  ni 
temperatures  on  a  permancr.t  file  (which  mi:  . 
used  as  input  to  subsequent  jobs)  and  use  u,  . 
results  to  perform  the  stress  analysis. 

a  2 

Run  the  heat  transfer  problem,  save  tin:  oh.-M/id 
temperatures  on  a  permanent  fiLe,  and  stop. 

=  3 

Run  the  heat  transfer  proLfl cm  and  stop. 

II- 6 
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C.  DATA  CARD  INPUT 

1.  Grid  Definition 

a.  Start  Cards 

(1)  First  Card  (A3) 
cc 

1-3  Enter  the  characters  nWAB 

(2)  Title  Card  (12A6) 
cc 

1-72  HED  Title  information  for  job  and  plot  (center  about  cc  36 
for  plot) 

b.  Control  Card  (One  card,  315) 


cc 

1-5 

NNP 

Number  of  nodal  points  <_ 

273 

6-10 

NEL 

Number  of  quadrilateral  elements  £ 

240 

11-15 

IPF 

Plot  flag  <  0,  plot  only 

=  0,  run  only 
>  0,  plot  and  run 

c. 

Node  Coordinate  Cards  (15,  2F10.3,  15,  2F10.3,  15) 

cc 

1-5 

N 

Node  Number 

6-15 

R(N) 

Radial  Coordinate 

(L) 

16-25 

Z(N) 

Axial  Coordinate 

(L) 

26-30 

NE 

Ending  Node  Point  Number 

31-40 

r(ne) 

Radial  Coordinate 

(L) 

41-50 

z-(ne) 

Axial  Coordinate 

(L) 

51-55  NI  Node  Number  Increment 
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If  the  ending  node  number  NE  is  not  zero  (or  blank),  then  nodes 
will  be  generated  in  equal  distance  increments  along  a  line  between  node  N  and  node  NE. 
The  first  generated  node  is  assigned  the  number  N  +  NI;  the  second  generated  node 
number  is  N  +  2NI,  etc.  Note  that  the  node  number  difference  (NE-N)  must  be  positive 
and  divisible  by  NI. 


If  is  omitted  it  will  be  assigned  a  value  of  *1*  automatically, 
d.  Element  Numbering  Cards  (8l5) 


cc 

1-5 

N  Element  Number 

6-io 

IX  (N,  1) 

r 

11-15 

IX  (N,2) 

Node  numbers  describing  the 
comer  points  of  the 

» 

16-20 

IX  (N,3) 

quadrilateral  ! 

i 

21-25 

IX  (N,4) 

s.  ^ 

26-30 

MN(N) 

Material  Number 

31-35 

NLG 

Number  of  elements  to  be  generated 

36-40 

NNI 

Node  Number  Increment 

If  the  number  of  elements  to  be  generated  (NLG) 

elements 

will  be  generated.  The  first  generated  element 

the  number  N  +  1;  the  second  generated  element  is  numbered  N  +  2,  etc.  The  node 
numbers  of  the  second  generated  element  are  found  by  adding  NNI  to  the  node  numbers  of 
the  first  generated  element,  etc.  If  NNI  is  omitted  it  will  be  assigned  a  value  of  *1". 


The  material  number  assigned  to  generated  elements  is  the 
same  as  that  for  element  N.  Material  numbers  are  not  necessary  if  the  program  is  to 
plot  the  grid  only  (lPF<  0). 


If  IPF  <  0;  data  input  ends  here. 
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Solution  Time 

and  Temperature  Information 

a. 

Control  Card  (One  card  215,  F10.0,  215,  F10.0). 

cc 

1-5 

NDT 

Total  number  of  time  increments  to  be  used 
in  the  solution  for  element  temperatures  and/or 
stresses. 

6-10 

NTR 

Number  of  regions  on  the  time  axis  having  the 
same  value  of  time  increment  ratio.  <  50 

11-20 

TMF 

i.  L 

Value  of  time  at  the  end  of  the  NDT —  increment,  (t) 

21-25 

NBCF 

Number  of  functions  describing  time  dependent  boundary 
conditions.  <  10 

26-30 

MEM 

Temperature  information  flag: 

=  0  Element  temperatures  are  to  be  calculated 

using  the  heat  transfer  analysis. 

=  1  Element  temperatures  are  to  be  read 

from  cards. 

=  2  Element  temperatures  are  to  be  read  from 

a  tape  created  on  a  previous  run. 

=  3  Element  temperatures  are  to  be  read 

from  cards,  and  this  distribution  is 
to  be  used  for  all  solution  time  'ointu. 

31-^0  TZ  Temperature  of  eveiy  element  at  time  zero  (*V) 

(stress  free  temperature). 
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Solution  Time  Points 

At  least  one  card  in  this  section 

First  Card  (3  (I5,2FLO.O) 


NTI(l) 


6-15  DT  (1) 


41-50 


Number  of  increments  for  which  the  time 
increment  ratio  C(l)  remains  constant  in 
region  1. 

Value  of  time  after  NTI(l)  time  increments 
in  region  1. 


16-25  C  (1)  Time  increment  ratio  in  region  1. 

26-30  NTI  (2)  Same  as  cc  1-5  for  region  2. 


31-40  DT  (2) 


66-75  C  (3) 


Value  of  time  after  NTI  (1)  HNT1  (2) 
time  increments  (T) 


C  (2)  Same  as  cc  16-25  for  region  2 


51-55  NTI  (3) 

56-65  DT  (3)  Region  3 


Second  Card  (3  (15,  2  F10.0))  (If  required) 


1-5  NTI  (4) 

6-15  DT  (4)  Region  4 

16-25  C  (4) 


Use  as  many  cards  In  this  section  as  are  required  to 
enter  NTR  groups  -£  (Nil  (I),  OT  «>.  C  W).  Note  that 

l  NTI  =■  NDT,  otherwise  an  error  message  will  be 
1=1  ^  issued  by  the  program. 


The  solution  time  points  in  region  r  are  found  using: 


'i+1  i 


-  t±  +  cr*  (t±  -  t^) 
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This  method  of  solution  time  point  input  has  the  effect  of  automatically 
producing  small  time  steps  at  the  beginning  of  the  region  while  continuously 
widening  the  time  steps  toward  the  end  of  the  region.  As  an  example,  if 
50  solution  time  points  are  specified  over  a  24  hour  period  with  C  *  1.1  the 
first  solution  time  point  will  occur  at  .0206  hrs  (1.24  minutes)  and  the  last 
solution  step  (t5Q  -  *-9)  will  be  2,19  hours  (132  minutes). 

It  is  important  to  choose  C  such  that  the  smallest  time  step  will 
not  be  so  small  as  to  produce  a  reduced  time  falling  off  the  mechanical 
properties  table  or  the  final  time  step  too  large.  The  following  equations 
may  be  used  to  find  these  values: 


C-l 

CN~1 


(Czll 

-  1 


where 


C  *  time  increment  ratio 
AT.^  =  Initial  solution  time  step 

AN*  total  number  of  solution  time  points 
ATN  *  final  time  step 

Tf  =  final  value  of  time  (X  *  0) 


(H-2) 


H-l 


These  two  equations  are  plotted  for  various  values  of  C  in  Figures 
and  H-2  respectively. 


H-ii 


)/(Time 
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L3 


Figure  H-2 
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c.  Boundary  Condition  Function  Cards 


Skip  this  section  if  NBCF  =  0 


(.1) 

Card  One  (215) 

cc 

1-5 

N  Number  assigned  to  this  : 

6-10 

NPTS(N)  Number  of  points  in  the  ■ 
this  function. 

(2) 

Card  Two  (SFIO.O) 

cc 

1-10 

TFN  (N,  l)Time  at  point  1,  t1  (T) 

11-20 

FN  (N,  l)  Value  of  the  function  at 

21-30 

TFN(N,  2)  Time  at  point  2,  tg 

31-40 

FN  (N,  2)  Value  of  the  function  at 

41-50 

51-60 

M  }  Point  3 

61-70 

71-80 

(3) 

Card  Three  (8F10.0)  (if  .required) 

cc 

1-10 

11-20 

(N,5)  1  Point  5 

FN  (N,5)  )  Plt  ' 

<  30 


(T) 


etc. 


Use  as  many  cards  2,  3>  etc.,  in  this  section  as  are  required 
to  enter  NPTS(N)  pairs  of  (TFN(N,  X),  FN(N,  I))  which  define  this  function  (Number  N). 
There  are  NBCF  sets  of  cards  1,  2,  3,  etc.  in  this  section.  Data  for  a  new  function 
begins  on  a  new  card  1. 


Figure  H-3  represents  a  function  which  might  be  used  to 
describe  the  pressure  transient  in  a  motors  The  last  point  in  the  table  must  have  a 
value  of  time  which  is  greater  than  or  equal  to  the  length  of  solution  period  (TMF); 
for  this  case  l6.0>  14.0  =  TMF  (Figure  3.l) . 
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The  table(s)  input  in  this  section  are  used  in  the  prescription 
of  time-dependent  boundary  conditions  such  as  environmental  or  nodal  point  temperatures, 
nodal  point  forces  and/or  displacements,  pressure  or  acceleration  loads,  etc.  One 
function  may  be  used  to  describe  multiple  types  of  boundary  conditions.  Values  of  the 
functions  at  times  other  than  TFN  (N,l)  are  calculated  within  the  program  using  linear 
interpolation. 

If  no  data  are  input  in  this  section  (NBCF  =  0),  then  all 
boundary  specifications  are  independent  of  time  and  boundary  values  are  assumed  to  be 
imposed  as  step  function  at  t  =  0*“ . 

d.  Printing  Control  (l6l5)  £  10  regions 


Card  1 

cc 

1-5 

6-10 

N0(l) 

NPR(l) 

Output  Interval 

Number  of  print  operations  at  this  interval 

11-15 

16-20 

N0(2) 

NPR(2) 

Output  interval 

Number  of  print  operations  at  this  interval 

71-75 

76-80 

•  •  ♦  • 

N0(8) 

npr(8) 

Output  Interval 

Number  of  print  operations  at  this  interval 

Card  2  (if 

required) 

CC 

I- 5 

6-10 

II- 15 
16-20 

N0(9) 

NPR(9) 

NO  (10) 

NPR  (10) 

Output  Interval 

Number  of  print  operations  at  tnis  interval 
Output  Interval 

Number  of  print  operations  at  this  interval 

This  information  is  used  to  control  the  amount^oT printed 

output  produced  by  the  program.  The  output  scheme  defined  in  this  section  is  used  in 
both  the  thermal  uud  stress  analysis  phases,  in  preparing  the  data  for  this  block, 
tl>e  following,  condition  mm  si.  be  observed: 

N0(l)  *  NPR(l)  +  N0(2)  *  NPR(2)  +  •••  —  NDT 
This  is  to  say  that  tae  printing  range  must  cover  the  solution  period. 
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Suppose  that  temperature  output  is  required  at  time  points  2,  4,  '( 
and  8  of  Figure  2.1,  then  one  card  would  be  input  in  this  section: 

Card  1:  2, 2, 3, 1,1,1  ;  [2  (2)  +  3  (l)  +  1  (l)  =  83. 

If  temperatures  at  all  8  points  is  requested,  then 

Card  1:  1,  8  ;  [l  (8)  =  8], 

e.  Average  Element  Temperature  Cards  (315,  F10.0) 

(Card  2-a  is  the  reference  for  values  of  NDT,  NTEM  and  TZ) . 


Skip  this  section  if  NTEM  =  0,  or  NTEM  =  2 


cc 

1-5 

NEIS 

Starting  element  number  for 

thfs  group. 

6-10 

NELG 

Number  of  elements  with  the 
as  element  number  NEIS. 

same  average  temperature 

11-15 

NELI 

Element  number  increment. 

16-25 

TAVG 

Aye rage  element  temperature 

for  this  group  of 

elements.  (°F) 

If  NIEM  =  1,  there  are  NDT  sets  of  element  temperatures 
that  must  be  defined  in  this  section  (one  set  for  each  of  the  NDT  time  points  contained 
in  the  solution  period).  All  element  temperatures  must  be  prescribed  at  a  given  time 
point  before  proceding  to  the  next  point.  If  every  element  has  a  different  temperature, 
then  HEL  cards  (with  cc  6-15  blank)  must  be  prepared  for  that  time  point.  It  is  possible 
to  generate  -element- 'temperatures  at  a  time  point  if  several  elements  arc  at  the  same 
temperature.  NEIS  elements  are  assigned  average  temperatures  TAVG.  The  number  assigned 
to  the  first  generated  element  is  NEIS  +  NELI;  the  second  is  NEIS  +  2  (NEI.I),  etc. 

Suppose  that  a  body  described  using  50  elements  (NEL  =  50) 
is  at  a  uniform  77°F,  when  t  <  0  (TZ  =  77.0),  and  at  the  end  of  the  fir  at  time  increment 
elements  1-25  drop  to  6Q°F  while  elements  26-50  reach  45°F.  Two  cards  are  required  to 
define  the  element  temperature  distribution  at  the  end  of  time  increment  one  (t  =  t-^) : 

Card  1:  1,  24,  1,  60.0 

Card  2:  26,  24,  1,  45.0 
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If  NTEM  =  3,  then  there  is  only  one  set  of  element  temperatures 
which  must  be  input  in  this  section.  This  temperature  distribution 
applied  to  all  solution  time  points;  the  body  starts  at  t  »  0”  at  a 
uniform  temperature  TZ  (at  which  the  body  is  assumed  to  be  "stress  free") , 
and  at  t  =  0+  the  ele-  "int  temperatures  assume  the  values  prescribed  in 
this  section  of  the  data  input  and  remain  constant  for  all  time, 
t  >  0. 

3.  Transient  Heat  Transfer  Solution  Data 


a.  Control  Card  (One  card;  415) 
cc 

1-5  NMAT  Number  of  materials  with  different  thermal  properties.  <  10 

6-10  NNBC  Number  of  nodal  point  boundary  conditions  (temperatures 

or  heat  fluxes) 

11-15  NCBC  Number  of  convection  boundary  conditions  <_  65 
16-20  JOB  Job  Control  Flag 

•*  0  :  run  heat  transfer  and  use  the  results  to  perform 
the  stress  analysis. 

=  1  :  run  heat  transfer,  save  element  temperatures  on 
a  permanent  file,  and  use  the  results  to 
perform  the  stress  analysis. 

=  2  <•  run  heat  transfer,  save  the  results  on  a 
permanent  file,  and  stop. 

=  3  :  run  heat  transfer  and  stop. 


21-25  KAT  KAT  =  0  Axisymmetric  Analysis 
KAT  i  0  Planar  Analysis 

b.  Material  Cards  (110,  6F10.0) 

cc 

I- 10  N  Material  number 

II- 20  XCOND(N)  Conductivity:  Krr,  (BtuMT)"1  (L)"2  (°F/L) 
21-30  YCOND(N)  Conductivity:  Kzz 
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31-40  XYC0ND(N) 
41-50  SPHT  (N) 
51-60  DENS (N) 
61-70  QX(N) 


Conductivity:  K 

J  rz 

Specific  Heat,  (Btu)  (FJ^C0?)"1 
Weight  density,  (F)  (L)  ^ 

Heat  generated  per  unit  volume,  (Btu)(T)  "'"(L) 


Use  one  card  for  each  different  material  number  assigned  in  the 
element  array  (cards  1-d) ;  NMAT  cards  must  be  prepared  in  this  section  - 
order  is  unimportant. 

If  the  thermal  conductivity  of  a  material  is  independent  of 
direction,  then 


K 


rr 


k 


rz 


K  *  k 
zz 

0 


The  heat  generated  per  unit  volume  is  assumed  to  be  constant 
with  time. 


c.  Nodal  Point  Boundary  Conditions  (215,  F10.0,  15) 


Skip  this  section  if  NNBC  *  0 


cc 

1-5 

N 

Node  Number 

6-10 

KODE(N) 

Boundary  condition  type 

f  *  0,  externally  supplied  heat  flux 
|*1,  prescribed  node  temperature 

11-20 

T(N) 

Boundary  value  amplitude 

f *  Heat  flux  (KOBE(N)  -  0),  (Btu) (T) 
\*  Temperature  (K0DE(N)  »  1),  (°F) 

21-25 

NFN(N) 

Function  number 

All  nodal  points  not  specified  in  this  section  are  assumed  to 
have  externally  supplied  heat  flux  of  zero  for  all  values  of  time. 

A  function  number  equal  to  zero  (or  blank)  means  that  the 
prescribed  boundary  condition  is  applied  at  time  zero  and  remains  constant 
for  ail  time,  t  >0. 
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The  functions  assigned  in  this  section  must  have  been  defined 
previously  in  Section  2-c.  For  time  varying  boundary  conditions,  the 
magnitude  of  the  boundary  value  at  some  time  t  is  found  by  selecting  the 
value  of  the  function  at  t  and  then  multiplying  this  value  times  the  boundary 
value  amplitude  (cc  11-20).  A  given  function  can  be  used  to  describe  any 
number  of  boundary  conditions. 

d.  Convection  Boundary  Conditions 


Skip  this 

section  if  NCBC  =  0 

(NCBC  Cards,  215 

,  2F10.0 ,  15) 

cc 

1-5 

I(N) 

Node  number  i 

6-10 

J(N) 

Node  number  j 

11-20 

H(N) 

Heat  transfer  coefficient  ,h:  (BTU)  (T)  ^(L)  ^  ("F)  * 

q  «  h  (T  -  Tq) 

21-30 

TE(N) 

Environmental  temperature  amplitude,  Tq 

31-35 

NFCV(N) 

Function  number 

If  the  environmental  temperature  T  is  time  dependent,  then  a 
non-zero  function  numer  must  be  specified  in  cc  3^-35.  Tq  will  be  mutliplied 
by  the  appropriate  value  of  the  function  of  time  t  in  order  to  establish  the 
value  of  environmental  temperature,  Tq. 

If  the  environment  does  not  change  temperature  with  time, 

then  NFCV(N)  =  0  and  T  *  T  ,  constant  for  t  >  0. 
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4.  Stress  Analysis  Information 


The  program  will  not  read  data  in  this  section  if 
NTEH  ■  0  (card  2a)  and  JOB  =  2  or  3  (card  3a) 


a.  Title  card  (20A4) 
cc 

1-8  HED  Any  alpha  numeric  information  (printed  with 

the  solution) 


b.  Control  Card  (515,  2(F10.0,  15),  315,  F10.0) 
cc 


1-5 

NMAT 

Number  of  materials 

£  4 

6-10 

NCMH 

* 

Number  of  elements  with  material  identification 
numbers  which  are  to  be  redefined. 

11-15 

NBCN 

Number  of  node  points  for  which  boundary 
cards  are  used 

1  60 

16-20 

NPC 

Number  of  pressure  cards 

1  55 

21-25 

NDMG 

Number  of  elements  for  damage  evaluation 

<  20 

26-35 

ANGV 

Angular  velocity  amplitude,  (R)  (T)  1 

36-40 

NFAV 

Function  number  for  ANGV 

41-50 

AZZ 

-2 

Axial  acceleration  amplitude,  (L)  (T) 

51-55 

NFAZ 

Function  number  for  AZZ 

56-60 

IDSF 

Pressure  boundary  condition  function  no. 
for  use  of  pressurization  shift  function 

data 

61-65 

IPSC 

Geometry  type  flag 

IPSC  ■  0  Axi'jymmetric  analysis 

■  1  Plane  strain  analysis 

■  2  Generalized  plane  strain  analysis 
*  3  Generalized  plane  stress  analysis 

66-70  NSM  Material  No.  (if  IPSC  =  2)  of  case  material 

71-80  SZV  f  Value  of  normal  stress  (if  IPSC  =  3) 

(F)  (L)-2 
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NMAT  is  the  total  number  of  materials  (viscoelastic  and 
elastic) . 


If  the  material  I.D.  numbers  assigned  to  the  elements  in  Section  1-d  are 
appropriate  for  both  the  heat  transfer  and  stress  analyses,  then  NCMN  is 
set  to  zero  (or  blank).  NBCN  is  a  count  of  nodes  at  which  force  and/or 
displacement  boundary  values  are  specified.  NPC  is  the  total  number  of 
element  sides  subjected  to  pressure  loads. 

IDSP  is  the  boundary  condition  function  number  which  describes  the  bulk 
pressure  as  a  function  of  time  when  an  isothermal  pressurization  case  is 
being  run.  For  a  thermal  analysis  leave  this  field  blank.  This  function 
will  be  used  to  find  the  pressurization  shift  function  values  in  the 
properties  determination  and  also  will  be  used  in  the  damage  calculations. 
When  IPSC  -  2  a  special  generalized  plane  strain  analysis  will  be  run. 

The  normal  strain  will  be  set  equal  to  the  average  thermal  strain  (a  AT  ) 
in  the  case  which  is  found  by  the  material  number  NSM  specified  in  cc  6§-70. 
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c.  Material  Properties 

(l)  Control  Card  (315,  4F10.0) 
ce 


1-5  K  Material  number 

6-10  NON(K)  Number  of  terms  in  the  Prony  series  representation  of  .1  36 
the  shear  relaxation  function. 


11-15  NSFP(K)  Number  of  points  in  the  shift  function  table.  1  16 

16-25  APO(K)  Equilibrium  shear  modulus,  (F)  (L)~2 
26-35  XK(K)  Bulk  modulus,  (f)(L)-2 

36-45  ALP(K)  Linear  coefficient  of  thermal  expansion,  (L)(L)~1(0F) 

Jl  p 

46-55  DENS(K)  Mass  density,  (F)(L)  (T) 


The  shear  relaxation  function  is  written  in  the  form: 

M  -At 

0  (t)  =  A  +  £  A,  e 
i=l 

where 

NON(K)  ~  M  for  the  K—  material 
APO(K)  =  Aq  for  the  k£"  material 


An  elastic  material  is  input  by  leaving  cc  6-15  blank  and  entering 
the  shear  and  bulk  moduli  in  cc  16-25  and  cc  26-35,  respectively.  A  non-zero 
value  of  density  is  required  for  the  calculation  of  body  forces  arising  from 
specified  values  of  spin  velocity  and/ or  axial  acceleration. 

(2)  Prony  Series  Coefficients  Card(s)  (8F10.0) 


Skip  this  section  if  NON(K)  =  0 


cc 

I- 10  AP(K,  1)  A-  for  material  K 

II- 20  BP(K,  1)  fi  ' 

21-30  AP(K,2)  Ap 
31-40  BP(K,2) 

41-50  AP(K,3) 

etc. 
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cc 

1-10 

FST(K,1) 

Temperature  at 

1st  point, 

(°F) 

11-20 

FS(K,1) 

Logl0  aTL 

21-30 

FST(K,2) 

Temperature  at 

2nd  point,  T2 

(°F) 

31-40 

FS(K,2) 

L°g10  aT 

41-50 

FST(K,3) 

<b 

etc. 

Use  as  many  cards  in  this  section  as  are  required  to  specify 
NSFP(K)  pairs  of  (T^,  log.^  aT  );  four  pairs  per  card.  If  IDSP  >  0  (card 

4b)  FST  (K,  1)  will  be  interpreted  as  pressures  and  FS  (K,  1)  are  Log1Q  ap. 

d.  Element  Material  Numbers 

Skip  this  section  if  NCMN  ■  0 


cc 

1-5 

N1 

(1) 

Element  number 

6-10 

N2 

(1) 

Material  number  assigned  to 

element  N1  (1) 

11-15 

N3 

(1) 

Number  of  elements  with  the 
as  element  N1  (1) 

same  material  number 

16-20 

N4 

(1) 

Element  number  increment 

21-25 

N1 

(2) 

Element  number 

26-30 

N2 

(3) 

Material  number  assigned  to 

element  N1  (2) 

31-35 

N3 

(2) 

Number  of  elements  with  the 
as  element  N1  (2) 

same  material  numbe*' 

etc, 

i 
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Use  as  many  cards  in  this  section  as  are  required  to  re-define  the 
material  numbers  of  NCMN  elements;  16  entries  per  card  are  possible. 

Suppose  that  a  50  element  problem  is  to  have  all  its  material 
numbers  changed  (NCMN  *  50) ,  and  all  even  numbered  elements  are  material  1 
while  all  odd  elements  are  material  2;  one  card  describing  these  changes 
would  read: 


[1,  2,  24,  2,  2,  1,  24,  2] 
e.  Node  Point  Boundary  Specification (s) 


At  least  one  card  in  this  section 


(15, 

2(15,  F10.0 

,  15)) 

cc 

1-5 

K 

Node  point  number  (NB(N)  *  K) 

6-10 

NFLR(N) 

Radial  boundary  condition  type 

{ 

=  0;  externally  applied  force 
*  1;  specified  displacement 

(F)  (R)"1 
(L) 

11-20 

BVR(N) 

Radial  boundary  value  amplitude 

21-25 

NFNR(N) 

Function  Number 

26-30 

NFLZ(N) 

Axial  boundary  condition  type 

1 

f 

*  0;  externally  applied  force 

*  1;  specified  displacement 

(F)  (R)"1 
(L) 

31-40 

BVZ(N) 

Axial  boundary  value  amplitude 

41-45 

NFNZ(N) 

Function  number 

46-55 

BUTH(N) 

Skew  boundary  angle  (degrees) 

A  total  of  NBCN  cards  must  be  prepared  in  this  section.  The  axial 
displacement  at  one  node  must  be  specified  as  a  minimum  requirement. 

Positive  boundary  values  are  in  the  direction  of  the  positive 
coordinate  axes. 


H-25 
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Zero  (or  blank)  function  numers  assigned  to  boundary  value  components 
implies  time  independence  (constant  for  t  >  0).  The  variation  of  a  boundary 
value  with  time  is  determined  by  multiplying  amplitude  times  the 
appropriate  value  of  the  corresponding  function.  All  nodes  not  specified 
in  this  section  are  assumed  to  have  no  externally  applied  loads  and  are  free 
to  displace  as  the  solution  dictates. 

The  shew  boundary  is  shown  in  the  figure  below.  If  BUTH(N)  +  0 
the  boundary  conditions  are  expressed  in  the  n-s  system. 


y.z 


f .  Pressure  Loads 


Skip  this 

section  if  NPC  -  0 

1 

J 

cc 

1-5 

1PC(N) 

Node  number  i 

6-10 

JPC(N) 

Node  number  j 

11-20 

PR(N) 

Pressure  amplitude 

(F>  (L)~2 

21-25 

NFNP(N) 

Function 

There  are  NPC  cards  in  this  section.  NFNP(N)  -  0  means  that 
pressure  is  applied  as  a  step  function  at  t  *  0*. 

Positive  pressure  acts  in  the  direction  shown  in  Figure  4.1. 


Aerojet  Solid  Propulsion  Company 
Report  1341-26F 


g.  Damage  Parameters 
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cc 

1-5 

LDMG(N) 

Element  no.  whose  damage  is  to  be  evaluated 

6-15 

STZR(N) 

ato 

16-25 

TZR(N) 

See  below 

0 

26-35 

SCR(N) 

There  are  NDMG  cards  in  this  section.  LDMG(N)  may  be  positive 
or  negative.  If  positive,  the  hoop  stress  will  be  used  in  damage  calculations; 
if  negative  the  maximum  principal  r-z  stress  will  be  used.  The  damage  is 
evaluated  using 


PEDt  - 


(o(t')  -  ocr  +  P(t'))B 
<»to  ■  °cr>B  ap  <p<t'» 


If  IDSP  »  0  (card  4b)  then  P(t')  *  0  for  all  times.  If 
IDSP  >  0  the  P(t’)  will  be  found  from  the  boundary  condition  function 
Indicated  by  IDSP. 


Data  input  ends  at  this  point 


Aerojet  Solid  Propulsion  Company 
Report  1341-26F 


Appendix  H 


IV.  PROGRAM  OUTPUT  DESCRIPTION 

Output  from  the  program  includes: 

1.  Input  geometry,  material  properties  and  solution  time  information 
in  a  self-explanatory  format. 

2.  Nodal  point  temperatures  at  the  specified  time-point  interval. 

3.  The  radial,  tangential,  axial  and  shear  stresses  and  strains 
and  the  average  element  temperature  at  the  specified  time-print  interval. 

4.  The  damage  rate  and  accumulated  damage  for  specific  elements 
at  each  time-print  interval. 

The  printed  program  output  consists  of  (1)  "echo"  reproduction  of  the 
data  input  cards  identified  in  self-explanatory  format,  and  (2)  results  of 
the  analysis  which  might  include  nodal  point  temperatures,  element  stresses/ 
strains,  node  displacements,  etc. 

A.  ELEMENT  NODE  NUM3ERINC 

It  should  be  noted  that  element  node  number  data  may  not  be 
listed  in  the  same  order  as  these  data  appear  on  the  input  card.  The  program 
logic  permutes  the  order  of  the  element  node  mmbers  so  that  the  largest  node 
number  is  always  last  (4th)  in  the  printed  list  for  each  element.  For  example, 
if  element  number  53  (shown  below) 


is  input  as  (1,  2,  21,  20),  the  program  will  print  the  data  as  (20,  1,  2,  21) 
so  that  "21"  is  last  while  the  original  counter  clockwise  order  is  preserved. 
Efficiency  is  gained  if  the  user  specifies  the  4th  node  as  the  largest  for 
every  element. 

The  reason  for  having  the  largest  node  as  last  one  in  the  sequence 
for  a  given  element  is  that  three  (3)  equations  are  assigned  to  this  node  as 
opposed  to  only  two  (2)  equations  per  node  for  the  others,  A  node  number  that 
is  not  the  largest  one  in  any  element  has  the  R  and  Z  displacement  components 
as  its  unknowns  (2  total).  All  other  nodes  have  R,  Z  displacement  components 
plus  the  "mean  pressure  variable"  H  as  unknowns  at  that  point  (3  total).  Thus, 
node  "21"  of  eleraent"53"  above  has  assigned  to  it  UR21>  VZ21*  a,‘fI  ,!53  as 
unknowns . 
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In  order 
!'H"  for  every  element, 


co  assign  equations  uniquely  to  the  unknown  values  of 
the  program  has  the  following  restriction: 


The  finite  element  quadrilateral  mesh  must  be  numbered 
so  that  any  node  number  is  maximum  in  only  one  (1) 
element. 


With  a  "layered"  numbering  scheme  (across  the  "narrow"  direction 
of  the  grid)  the  restriction  is  complied  with  automatically: 


The  "heavy"  nodes  (•)  shown  above  have  three  (3)  unknowns 
associated  with  them. 


B.  SIGN  CONVENTION 

All  positive  vector  quantities  (forces/displacements)  are  in 
the  directions  of  the  positive  coordinate  axes  providing  the  grid  is  positioned 
in  the  "positive"  R-Z  quadrant: 


y.z 


>  x,R 


Normal  stress  (strains)  are  tensile  when  positive  and  the  shear 
stress  (strain)  sign  convention  is  shown  below: 


x,R 


Aerojet  Solid  Propulsion  Company 
Report  1341-26F 


Appendix  H 

dt* 

If  IDSP  -  0  (Card  4b)  then  P(t')  *  0  for  all  times.  If  IDSP  >  0  the  P(t') 
will  be  found  from  the  boundary  condition  function  indicated  by  IDSP. 


PED.  »  «*- 
C  St 


t  (o(t»)  -  ocr  +  P(t’))B 


to 


-  acr)B  a  (P(f)> 


The  value  of  P(t')  is  that  portion  of  the  inner-bore  firing 
pressure,  P.(t),  which  is  transmitted  to  the  point  in  the  grain  where  the 
principal  stress,  o(t'),  is  evaluated.  The  relation  of  P(t’)  to  P. (t') 
in  infinite-length  cylinders  is  well  known (  )  (See  Section  IV  of  this 

report).  For  finite  length  cylinders  the  ratio  of  P(t')  to  P  (t  )  must 
be  obtained  from  the  stress  analysis.  *• 

For  metal  cases  P(t')  seldom  differs  from  P  (t1)  by  more  thai;  5Z. 
For  most  firing  problems  this  difference  is  not  significant  and  P(t') 
can  be  approximated  by  P  (t1). 


Data  input  ends  at  this  point 
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A  positive  spin  velocity  produces  hoop  tension,  and  positive 
axial  acceleration  causes  body  forces  to  be  applied  on  the  body  acting  in 
the  (-Z)  direction: 


Body  forces  resulting  from  a 
positive  Z-acceleration 


x,R 


V.  EXAMPLE  PROBLEMS 

A.  TRANSIENT  HEAT  CONDUCTION  IN  A  LONG  CYLINDER 

The  purpose  of  this  example  is  to  illustrate  the  use  of  the  program 
in  solving  heat  conduction  problems.  A  long,  hollow  cylinder  constructed  from+ 
a  single  material. is  initially  at  a  uniform  temperature  of  0°F;  then,  at  t  *  0 
the  outer  surface  of  the  cylinder  is  instantaneously  heated  to  1#F  along  its 
entire  length.  The  ends  of  the  cylinder  are  insulated  against  axial  heat  flow 
so  that  at  any  axial  station  the  heat  flow  is  purely  radial. 

Figure  H-4  shows  a  one  (1)  inch  slice  of  the  cylinder  whose  inner 
and  outer  radii  are  1.00"  and  2.00",  respectively.  The  slice  has  been  modeled 
with  a  mesh  composed  of  ten  (10)  equal  sized  elements;  the  "Z"  axis  is  the  centei — 
line  of  the  cylinder.  At  time  zero  the  temperature  of  nodes  11  and  22  is  changed 
from  0®F  (TZ  =  0.0)  to  1°F,  and  as  time  proceeds  the  interior  of  the  cylinder 
begins  to  warm  up;  thermal  equilibrium  is  reached  when  the  entire  cylinder  is 
at  a  uniform  1°F.  Nodes  1-10  and  12-21  are  insulated  in  the  sense  that  no 
externally  supplied  heat  enters  the  body  at  these  points.  The  program  assumes 
that  all  nodes  not  specifically  included  as  boundary  condition  nodes  are  insulated 
against  externally  supplied  heat  flow;  i.e.,  for  all  non-bcundary  nodes,  the 
amount  of  heat  entering  a  node  must  balance  the  amount  of  heat  leaving  that 
node  in  a  unit  of  time. 

The  solution  span  Is  sub-divided  into  fifty  (50)  equal  time 
increments  of  0.01  hours  each  so  that.the  time  at  the  end  of  solution  is 
0,050  hours  (TMF  -  0.05),  Figure  H-5. 
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R  *  1.0" 


T 

7.  =  2.0" 

~T 

Z  =  1.0" 


a.  Mesh  For  Example  Problem  A 


0.00 

^ - 

-  50  intervals  at  0.001  hours 

a.  Solution  Time  Point  Array 


TMF  *  0.05 

__4 - 


Time 


0  10 

j - + - 

— *|  10  @  1  \+ -  20  @  2 


50 

-4 - 

*1 


b.  Solution  Output  Intervals 
Figure  H-5:  Solution  Time  Point  and  Output  Schedules 
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Figure  H-5 
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The  output  printing  schedule  is  set  up  so  that  results  are 
printed  at  every  solution  time  point  for  the  first  ten  (10)  increments  and  then 
at  every  other  solution  for  the  remaining  forty  (40)  increments.  Figure  H-5. 

10  print  operations  at  an  interval  of  1  »  10 
20  print  operations  at  an  interval  of  2  *  40 

50  increments 


The  thermal  properties  of  the  material  are  summarized  as  follows: 

Krj.  -  0.20  (Btu)(hr.)“1(in.)"2(°F/in.)"1 

K  «  0.20 
zz 

K  *  0. 
rz 

Specific  Heat  *  0.20  (Btu)  (lb.)'1  (0F)_1 
Weight  Density  *  0.20  (lb.)  (in.)  ^ 

Heat  generated  per  unit  volume  per  unit  time  =  0.* 


The  data  cards  for  this  job  are  shown  in  Table  H-3. 

Figure  H-6  is  a  plot  of  temperature  versus  time  for  three  (3) 
points  in  the  cylinder:  (1)  outer  surfac"  (nodes  11  or  22);  (2)  mid-radius 
(nodes  6  or  17);  and,  (3)  inner-radius  (nodes  1  or  12).  For  long  times, 
all  nodes  approach  1°F  in  the  limit. 

A  value  of  3  was  assigned  to  the  control  variable  '‘JOB" 
causing  a  termination  of  execution  after  the  heat  trarsfer  solution. 

If  the  results  were  to  be  saved**,  then  the  node  temperatures  would  be 
averaged  for  each  element  before  saving  the  results. 

B.  ELASTIC  RING  SUBJECTED  TO  AXIAL  TENSION 

The  purpose  of  this  example  is  to  illustrate  the  use  of  the 
program  in  solving  elastic  problems.  A  hollow,  short  cylinder  is  subjected 
to  an  axial  tension  of  6000  psi  on  one  surface  and  restrained  (without  radial 
shear)  against  axial  displacement  on  the  opposite  end  (Figure  H-7).  The 
problem  has  been  modeled  with  four  (4)  quadrilateral  elements  as  snown  in 
Figure  H-8.  The  applied  stress  has  been  converted  to  equivalent  concentrated 
loads  of  4000,  12000,  and  8000  lbs. /radian  applied  at  nodes  7,  8  and  9, 
respectively. 

Output  from  the  computer  program  is  shown  in  Figure  H-9  a,  b  and  c. 

:  For  purposes  of  data  preparation,  unit  heat  generation  is  treated  as  a 
physical  property  which  can  vary  from  material  to  material. 

**J0B  =  1  or  2  x 
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FORTRAN  STATEMENT 


Temperature 


Figure  H-6 


EXAMPLE 
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11:  AXIAL  LOADING  OF  AN  ELASTIC  CYLINDER 


£ 


6000  psi 


6000  psi 
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GRID  FOR  EXAMPLE  B 


R-1.0  R=2.0  R=3  ,C 


Figure  H-8 
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The  solution  time  span  which  has  no  meaning  when  constant  loads  are  applied  to 
elastic  media  consists  of  one  (1)  time  increment  arbitrarily  selected  as  being 
60  seconds  long;  one  print  operation  is  performed  at  the  end  of  the  first  (and 
only)  time  increment.  The  average  temperature  of  all  elements  is  read  in  as  0°K. 
The  pertinent  mechanical  properties  are  the  shear  (G  =  6(10)6  psi)  and  bulk 
(K  *»  10(10)6  psi)  moduli  of  the  material.  There  were  eight  (8)  of  the  nine  (9) 
nodes  desig:..  -?d  as  boundary  nodes;  nodes  A  and  6  are  unloaded  and  ui  **es trained 
and  do  not  have  to  appear  in  the  statement  of  boundary  conditions.  Nodes  A  and  6 
were  inadvertently  retained  from  a  previous  data  deck. 


The  results  consist  of  R  and  Z  displacements  of  each  node,  average 
element  stresses  and  strains  and  the  average  element  temperatures  all  quoted  at 
t\  =  60  seconds.  The  applied  axial  stress  is  recovered  exactly  as  6000  psi  (SIGZ) 
con;  tant  in  each  element. 


C.  PLANE-STRAIN  V1SC0EIASTIC  CYLINDER,  EXAMPLE  PROBLEM 

The  problem  is  a  long  cylindrical  bore  propellant  grain  (1.60  I.D., 

8.0  O.D. )  in  a  .060  steel  case.  The  element  geometry  is  shown  in  Figure  H-10. 

The  propellant  relaxation  and  shift  function  curves  are  shown  in  Figu; es  H-ll  and 
H«?12.  Sixteen  decade  reduced  time  points  were  used  to  perform  a  prony  series 
curve  fit  to  the  relaxation  function  and  the  coefficients  are  shown  in  Table  i!— A . 

The  motor,  initially  stress  free  at  135°F,  was  subjected  to  the 
thermal  environment,  through  a  heat  transfer  coefficient  at  the  case,  shown  in 
Figure  H-13.  All  nodes  in  the  system  were  fixed  in  the  axial  direction  to 
simulate  a  plane  strain  condition  and  the  linear  cumulative  damage  was  calculated 
for  two  elements  in  the  system,  one  at  bore  and  one  next  to  the  case  -  grain  bond. 

A  listing  of  the  data  input  cards  required  to  describe  this  problem 
is  shown  in  Table  H-5.  Line  2  is  the  title;  line  3  is  the  control  card;  lines 
A-9  describe  the  node  point  coordinates;  lines  10,  11  establish  the  element  I.D.; 
line  12  is  the  solution  time  point  and  reference  temperature  control  and  lines 
13-29  describe  the  solution  time  points.  Lines  31-36  describe,  the  free  stream 
temperature  function.  Line  37  describes  the  output  print  interval.  Lines  39 
and  A0  are  the  thermal  properties'  for  the  grain  and  case.  Line  A1  describes  the 
convection  boundary  condition.  Line  A2  is  the  stress  analysis  title  card  and 
A3  is  the  stress  control  card.  Lines  A5-A3  are  the  relaxation  series  coefficients 
Lines  A9-51  are  the  shift  function  points.  Line  52  describes  the  case  properties. 
Lines  53-92  are  the  displacement  boundary  conditions.  The  last  two  lines  give 
the  damage  parameters. 

Selected  portions  of  the  output  from  this  analysis  are  shown  in 

Table  H-6. 


H-A2 
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Reduced  Time,  t/a  ,  min 


5 


- 
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Propellant  Relaxation  Function 

log(t/a_) 

(mih) 

E 

(pit) 

(“r1) 

Ai 

Ai/3 

i 

100 

33.333 

0 

-9 

44000 

3  x  1010 

25836, 

8612 

1 

-8 

29000 

3  x  109 

12087 

^  '4029 

2 

-7 

14000 

3  x  108 

6919.2 

2306.4 

3 

-6 

3300 

3  x  lO7 

3915.5 

1305.2 

4 

-5 

5000 

3  x  106 

2435.2 

8li,40 

5 

-4 

3000 

3  x  105 

1341*9 

447.3 

6 

-3 

1800 

,3  x  104 

1033.9 

344.63 

7 

-2 

1000 

i  x  io3 

457.98 

152.66 

8 

-1 

620 

3  k  102 

265.80 

88.60 

9 

0 

410 

3;  x  101 

137.74 

45.91 

10 

1 

320 

3  x  106 

.16481 

.0546 

n 

2 

280 

3  x  10"1 

113.47 

37.82 

12 

3 

210 

3  x  10“2 

r-2.9029 

-.9676 

13 

4 

190 

3  x  10t3 

60.2506 

20.0835 

14 

5 

150 

3  x  10"4 

4.26975 

1.42325 

15 

6 

130 

3  x,  ;10"5 

49.4143 

1,64714 

16 

IS  -B  t 

E  -  A 

r  n 

+  l  A,  e  1 

a  v 

1-1  1 

H-46 

Table  H-4 
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Figure  H-.1/J 
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Appendix  1 

NON-LINEAR  ANALYSIS  BASED  ON  PROPELLANT  DILATATION 

In  general,  the  response  of  grains  to  various  loading,  conditions 
is imainly  determined  by  the  dilatational  behavior  of  the  propellant,  hence, 
it  would  appear  that  accounting  for  the  non-lihear  dilatational  response 
of  propellant  is  of' primary  importance. 

Although  insufficient  experimental  data  is  available  to  permit  a 
comprehensive  characterization  of  the  non-linear  dilatational  effects ,  it 
'was  felt  that  an  evaluation  could'  he  obtained  by  utilizing  a-  description 
which  qualitatively  accounts  for  the  behavior  reported  for  a  limited  class 
of  stress  states  in  References  i-1  and  1-2,  and  which  predicts  reasonable 
results  for  other  stress  states.  Because  of  the  tentative  nature  of  the 
characterization,  it  was  deemed  desirable  to  approximate  the  rion-linear 
behavior  by  an  equation  which  could  be  incorporated  into  the  existing 
analysis,  with  a  minimum  amount  of  effort. 

With  the  above  objectives  in  mind  the  following  approximations 
were  made:  £1)  The  thermo-rheologically  simple  linear  viscoelastic 
distortions!  stress-strain  law  and  the  temperature  induced  volume  change 
relationships  Would  remain  unchanged;  (2)  the  stress  induced  dilatational 
response  would  be  approximated  as  elastic.  It  was  felt  that  the  visco¬ 
elastic  dilatation  (compare  Figures  11  and  12  of  Reference  It2) ,.  is  of 
secondary  importance  as  compared  to  the  basic  non-linear  effects;. 
Additionally,  the  lack  of  experimental  evidence  concerning  dilatation 
during  unloading  precludes  its  -consideration  at  this  time.  (3)  The 
stress-dilatation  relationship  would  be  considered  to  be  temperature 
Independent  (it  was  felt  that  the  actual  temperature  dependence  was  of 
secondary  importance).  (4)  The  characterization  would  be  limited  to 
relatively  small  strains.  (5)  The  elastic  relationship  between  the 
dilatation  and  the  stress  state  could  be  expressed' as : 

(6  -  3a AT)  -  f(g(HH))  =  £(«r  )  (1-1) 

where  "f"  is  a  non-linear  algebraic  function  of  the  quantity  "g",  g  (M^) 
is  in  turn  a  non-linear  function  of  the  stress  invariants..  The  stress 
Invariants!  are  defined  by  the  equations 
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X  Vol.  %  *  %  filled  polyurethane  rubber 

a  ■  hydrostatic  pressure,  psi 
e  =  strain  rate<,  iii./in. 
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Figure  1-2.  The  form  of  the  function  g  (  >6^)  was  detemined  by  considering 
the  data  presented  in  Figure  1  of  Reference  (1-2).  A  number  of  different 
functional  forms  were  investigated  from  among  which  the  following  form 
was  selected. 


g(  *0^)  -•»*  -  Kj_  ^  -fc  K2  *0^  +  iC3  &  3  (1-6) 

The  values  of  K^,  and  K3  were  selected  by  noting  that  Equation  (1-1) 
states  that  for -a  given  amount  of  dilatation  the  value  of  "g"  should 
be  a  constant.  Thus  for  the  different  stress  states  (corresponding 
to  different  values  of  hydrostatic  pressure)  given  in  Figure  1  of 
Reference  (1-2)  an  attempt  was  made  to  select  the  function  "g"  so 
that  for  a  given  value  of  0  all  states  would  yield  a  common  value  of  ''gn. 
The  following  values  were  found  for  the  parameters 

Kx  «  4.9  x  l(f3 

K2  -  2.3  x  10"5 

K  -  5,7  x  10~6 
3 

Using  these  parameters  the  following  table  was  constructed 


.  i  .  -  :  “  .'*■ 

’  2 

3:‘ 

4 

5 

Hydrostatic 

Pressure, 

0  - 

.005 

0  =  .01 

psi 

8. 

0/3 

sT 

"m 

15 

66 

84 

68 

65 

66 

-5 

99 

14 

165 

lOl 

-150 

125 

-115 

If  the  dilatational  response  could  be  described  by  a  linear  elastic 
law  all  the  entires  in  Column.  3  (and  in.  Column  5)  would  need  to  be 
approximately  equal,  the  extreme  grossness  of  this  approximation  is 
evident.  If  the  non-linear  function:  "g"  is  to  represent  the  experimental 
data. all  the  entires,  in  Column  2  (and  in  Column  4)  would  need,  to  be 
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approximately  equal,  while  this  is  not  true,  the  very  substantial 
improvement  of  the  non-linear  representation  as  compared  to  the  linear 
representation  is  evident. 

While  no  experimental  information  was  available  for  stress  states 
other  than  the  uniaxial  test  with  superimposed  hydrostatic  pressure,  it 
is  extremely  important  that  the  proposed  model  predicts  reasonable  results 


for  other  stress  states  as  these  stress  states  may  occur  in  a  given 
grain  problem.  Utilizing  the  proposed  non-linear  model  dilatation  was 
predicted  for  the  following  stress  states: 

(a) 

Uniaxial  stress  with  superimposed  hydrostatic 
pressure  (Figure  1-3) 

(b> 

Hydrostatic  pressure  (Figure  1-4) 

(c) 

Biaxial  stress  with  superimposed  hydrostatic  pressure 
(figure  1-5) 

(d) 

Simple  shear  with  superimposed  hydrostatic  pressure 
(Figure  1-6) 

The  straight  lines  in  Figure  1-3  indicate  the  predictions  of  the 
linear  model,  the  curved  lines  the  predictions  of  the  non-linear  model. 

The  predictions  of  the  non-linear  model  qualitatively  appear  to  be  reasonable. 

The  Incorporation  of  the  non-linear  dilatational  relationship 
(Equation  (x~l))  into  'jhe  existing  thermoviscoelastic  analyses  presented 
some  difficulties.  The  existing  thermoviscoelastic  analysis  employs 
the  following  incremental  relationship  for  dilatation  (for  time  increment 
N) 


(A9N  -  3a A  Tn) 


*N 


+  X 


N 


fcN-l  - :t  -  fcN 


(1-7) 


In  the  previous  linear  analysis  was  considered  to  be  a  constant  (the 
elastic  bulk  modulus)  and  xN  “  0...  The  utilization  of  Equation  (t*7)  as 
ah  approximation  to  Equation  (1-1)  was  accomplished  as  follows: 


Having  obtained  the  solution  for  time  t^^  consider  the  utilization 
of  Equation  (1-7)  for  -  t^.  Symbolically  acinoting  the  stress  state 
at  time  t„  ,  as  t. .  ana,  the  incremental  change  in  stress  state  for 
N"1  '  1JN-1 


i  as  . .  . 

N-l  ljfj-i 


The  following  two  associated  stress  states  are  defined 


1-7 
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ATIOXAL  PREDICTIONS  FOR  BIAXIAL  STRESS  KITH  SUPERIMPOSED  HYDROSTATIC  PRESSURE 
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Tija 


(1-8) 


4* 


N-l 


(1-9), 


Utilizing, Equations  (1-1),  (1-2),  (1-8),  and  (l-9):y  the  following  quantities 

are  calculated  9  ,.  0^, >0^  and  >0^,  The  following  approximate  relationship 

is  written  for  values  of  0  between  0  and  0, 

•a  b 


0 


•0«  0,  -  *01  0 
a  b  b  a 

ra 


*0*  -'»9^ 


(1-10) 


Assuming  that  the  absolute  value  of  the  change  in  stress  during  the 


increment  At^  will  be  of  the  saane  order  of  magnitude  as  |At 


ij 


N-l 


(i,e. ,  using  Equation  (I  -10)  and  writing  0^  “  ®n-1  +  an<* 
^  +  A-  *0^  yields 


0N~1  +  A6N 


i^0,  - 

a  b  ... 


KR  -  *0^ 


0-0. 
a  b 

■A  i  |A  ( 

iv^  -  rv|^ 


$0^  +  ,A  *0j,)  (1-11) 


Comparing, the  above  expression,  to  Equation. (1-7)  yields 

0  -  9, 
if  =*  _a  Jl 

^  5*0«  -  10^ 


Vri9 

9a  ^b 


9  -  0. 

a  a.  b 

(0t  - 
a:  b 


^N-l  "  8N-1 


Because  of  the  secant  approximation  to  the  non-linear  function,  in  general, 
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BAjIC  cumulative  damage  equations 


The  LCD  relation  merely-  states  the  manner  in  which  damage  can  -be  added. 

To  apply  the  relation  in  a  three-dimensional  stress  field  requires  a  time- 
dependent,  stress  failure  criterion;  An  empirical  approach  showed  that  the 
applicable  relation  is  a  -time-dependent,  maximum  principal  stress,  failure 
criterion  (MPS). 

Combining  the  LCD  and  MPS  criteria  leads  to  integral  relations  which 
take  full  account  of  the  past  stress-time-temperature  history  at  a  point 
in  the  grain. 

1 .  Linear  Cumulative  Damage  Relations 

The  very  simple  relation,  applicable  to  solid  propellant  failures,, 
is  defined  in  terms  of  tests  made  under  a  constantly  imposed  "true"  stress,  ot 
For  a -number,  N,  of  discrete  stress  levels  the  accumulated  damage  fraction, 

ED,  is  given  by  the  following  linear  relation: 


N. 

E 

i«l 


(J-l) 


where  ED  is  the  cumulative  damage 

P(h)  is  a  statistical  distribution-  parameter  and  relates 

the  n™  test  specimen  in  the  distribution  to  the  mean 
•of  the  population. 

■  th 

At^  is  the  increment  of  t-ime  the  specimen  is  exposed  to  the  1 
"true"  stress  level 

tfi  is  the  mean  time-.to-failure  for  the  population  of  specimens 
if  the  specimens  saw  only  the  ic-  "true"  stress  level. 


The 

damage  required 


accumulated  damage,  ED,  thus  gives  the  fraction  of  that 
to  fail  the  specimen.  Thus,  by  definition  ED  =  1  at  failure. 


The  parameter,  P(n),  is  a  highly  useful  term  in  that  it  provides 
the  focal  point  for  all  -the  statistical  studies  pertinent  to  the 
cumulative  damage  relations.  In  its  simplest  -forra--P(h-).  defines  the  position 
of  the  indiviaul  failure  with  respect  to  the  mean  of  the  distribution  of 
these  failures.  This  is  seen  on  considering  specimens  held- -under  a  single 
stress ,  where  P(ri)  equals  the  ratio  of  the  time- to-faiTure,  tf  (n),(  for  the 
ritn  individual  divided  ‘by  the  mean  timertorf allure  for  the  entire  -population, 


J-l 


A  — 


r 
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Thus,  we  have 


P-(n)  =  tf.(n)/t^. 


U-2) 


The  cumulative  damage  relation  merely  preserves  this  concept 
of  P(n)  for  conditions  of  complicated  loading  histories;, 

In  Reference  Jrl,  it  was  shown  that  P(r. )  is  independent  of  the 
stress  level.  This- justifies  placing  the  parameter  outside  the  summation 
in  Equation  (J-l) .  : 


Statistical  evaluations  of  P(h)  obtained  from  solid  propellant 
failure  testing,  show  the  parameter  to  follow  approximately  a  log  normal 
distribution. 


2 .  Time-Dependent  Maximum  Principal  Stress  Failure  Criterion 


The  most  general  form  of  the  MPS  failure  criterion  for  solid 
propellants  appears  to  be  a  modified  power-law  relation.  However,  a 
simple  power-law  approximation  to  this  criterion  covers  all  of  the  practical 
motor  problems  that  we  have  met  to  date.  This-. power-law  relation  has  the 
following  form: 


cf  -  Vt 


l  G . 

-  0 

t 

cr 

o' 

O' 

V  t. 

-  -cr 

-B 


(J-3)' 


where!  is  the  "true"  stress  applied  to  the  specimen. 

t^  is  the  mean  tirae-tor-failure  of  specimens  sheld  under  the  constant 
true  stress,.  ofc;. 

t  is  the  unit  value  of  the  time  for  whatever  units  are  used 
in  measuring  t^. 

atQ  is  the  true  stress  required  to  fail  the  specimen  in  the 

time  t  . 
o 

a  is  the  critical  true  stress  below  which  no  failures  are 
observed, 

a,p  is  the  time-temperature  shift  relation: 


j-2 
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It  has  been  found  to  apply  to^  solid  propellant  behavior 
with  correlation  Coefficients  usually  exceeding  0.93.  = — - 

The  limiting  stress,  a  ,  is  difficult  to  evaluate  because  it 
is  usually  not  large  compared"  with  the  relatively  large  variability-  typical 
of  solid  propellants.  It  appears  to  be  negligibly  small  (except  for  the 
effects  of  pressure)  for  all  the  propellants  tested  to  date. 

The  cumulative  damage  relation  is  generalized  as  follows. 

Firs*-,  Equation  (J-3)  is  combined  with  Equation  (J-l)  to  obtain  the 
cumulative  damage  relation  in  terms  of  discrete  stress  levels,  loading  time 
intervals,  and  test  temperatures. 


(J-4) 


where  a  .  is  the  i  stress  level. 

'  -  ti 


For  continuous',  changes  in  the  stress  and  in  the  temperature. 
Equation'  (J-4)‘  becomes : 


- - - 1 - j. — /fc 

P(n)  .  ccr)B  0 


(o  -  a  _  ) 
t  cr 

a,r(t) 


fl- 5-) 


where  ot>.  the  true  stress,  is  now  a  function  of  time, 

a,j,(t)  is  the  time-temperature  shift  relation  with  the  temperature 
expressed  as  a  function  of  time. 


Equation  (J-5)  represents  a  general  form  of  the  linear  cumulative 
damage  relation  for  solid  propellants.  This  equation  permits  the  summation 
of  damage  for  any  type  of  thermal  or  mechanical  loading  history,  provided 
the  stresses,  times  and  temperatures  are  known.  Also,  Equation  (J-5) 
tan  be  applied  without  change  to  three  dimensional  stress  problems; 


J-3 
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STUDY  OF  PROPELLANT  FAILURE  UNDER'  PRESSURE, 


A.  FAILURE  MECHANISMS; 

The  basic  mechanisms'  of  propellant  failure  have  been  extensively  studied 
undercurrent  and  oast  research  programs  (*'■-■1,  K-2),  Failure  is  considered 
to  be;  a  three  step  process  beginning  with  vacuole  formation.  The  vacuoles 
•may  initiate  in  che  polymer  or  at  the  binder-filler  interface,.  The 
associated  volume  change  may  be  so  negligible  that  the  process  is  termed  a 
"Mullin’ s  Effect".  The  second  step  is  the  extension  of  the  vacuoles  until 
they  overlap  producing  a  thin  membrane  of  polymer ^between  them.  The  third 
and  final  step  is  the  production  of  a  critical,  biaxial,  stress  field  in 
one  or  more  of  the  thin  membranes.  The  resulting  tears  propagate  vary 
slowly  relative  to  the  speed  of  sound  in  the  propellant.  This  failure 
mechanism  is  the  basis  of  the  considerations  made  later 

The  maximum- principal  stress  failure  criterion,  Equation  (J~3), 
contains  a  simple  stress  difference  term,  ot  -  a ^ .  Some  feel  that 
this  indicates  a  maximum  shear  stress  criterion  to  be  operative.  In 
biaxial  tension  this  is  expressed  by 

i (critical)  =  (o^  -  (Krl) 

where  r(critical)  Ir.  the  maximum  shear  strength  of  the  material 

0^  and  o ^  are  principal  stresses- 

Considering  the  stress  difference  terms  only,,  for  both  failure  criteria, 
leads  to  the  tabulation  in  Table  K-l.  It  is  assumed  that  a  maximum 
principal  stress,  0.,  is  imposed  in  tension.  The  remaining  stresses  are 
tabulated. 

Table.  K-l  shows  that  the'  two  criteria  give  identical  results  for 
simple  uniaxial  and’  strip-biaxial  tensile  data  from  which  the  MPS  criterion 
was  derived.  However,  the  two  criteria  differ  greatly/ when  applied  to  the 
<triaxial„  tensile,  stress  field  of  the  :poker  chip  specimen.  In  this  case, 
the,  stress  difference  for  ;the  maximum  shear  stress  criterion  is  dependent 
upon  1-k,  which  Is  very  small,  indicating  that  very  large  axial  stresses, 

0,  would  be  required  to;  fail  the  specimen.  On  the  other  hand,  the  MPS 
criterion  gives  a  relatively  large  value  for  the  -stress  difference  term, 

•0.  -  0  ..  This  difference,  like  those  for  the  uniaxial  and  biaxial 
tensile- data,  equals  the  axial  stress  on  the  specimen.  Tests  on  poker  chip 
specimens  showed*  the  MPS  criterion  to  hold, while  the  maximum  shear 
stress,  criterion  is  grossly  in, error. 


K-l 
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TABLE-  K-1 

STRESS  DIFFERENCES  AT  ATMOSPHERIC  PRESSURE 


MFS  Criterion 

'Maximum  Shear 

Stress  Criterion 

Test  Mode 

ft 

o  * ** ***  '  a  .  r  o 

cr  t  cr 

o^  a  2  o 

1.  ~  °2 

Uniaxial  Tensile 
Test 

a 

6  a 

o  0 

d 

Strip-Biaxial 
Tensile  Test 

a 

0  a 

a  0 

o 

Poker  Chip 

Tensile  Test 

a 

0  a 

a  ko 

(l-k)o 

*  qcr  is  set  equal  to  zero,  which  appears  to  be  thecaoe  for  many 
composite  propellants. 

**  The  maximum  tensile,  stress^  at  the  center  .of  .the  specimen. 

***  'k  Is  a  proportionality  factor  which,  in  this  case,  is  close  to  one, 
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This  result  is  not  surprising  in  view  of  the  fialure  mechanisms 
summarized  above.  Tensile  failures  in  uniaxial,  biaxial,  and  poker  chip 
specimens  follow  the  same  mechanism.  Once  the* vacuoles  are  produced' all  three, 
test  modes  become  equivalent,'  in  that  failure  is  produced  by  extension  of 
a  foam-like  material  until  the  biaxial  stress  in  the  membranes  exceeds 
the  critical  value  and  tearing  starts. 

B.  EVALUATING  THE  TIME-PRESSURE  SHIFT  FACTOR 

The  time-pressure  shift  factor  for  propellants  is  best  understood  on 
referring  to  the  foam-like  failure  behavipr  mentioned  above ,  Before  the 
vacuoles  are  formed,  or  when  they,  are  still  quite  small,  a  equals  one 
(up  toahout. 1500  psi  pressure),  This  was  shown  for ^tensile  moduli  by 
Wiegand,  '  Hazeltori'*'""^?  and  Lim  and  Tshoegel^ 
the  a  values  are  those  for  the  solid"  binder  which 
pressBre  changes  to  effect  significant  changes  In 
On  the  other  hand,  propellant  failure  data  are  sen 
pressure'  changes.  Dilatation  constitutes  the  only  physical  difference  in 
the  propellant  between  the  points  where  initial  tensile  moduli  and  failure 
data  are  obtained.  At  if allure,  this; dilatation  is  extensive  and  can  be 
partially  .suppressed  by  superimposed  hydrostatic  pressures.  As  shown 
by  Farris'  and  Lim  and  Tschoegel-™~5\  the  dilatation,  at  failure  is 
positive",  even,  at  the  higher  pressure  levels. 

Taken  together,  these  observations  indicate  that  a  provides  a  measure 
of  the  suppression  of  vacuole  forinaticin.  But,  p  is  Blso  a  measure  of 
the  same  affect.  That  is,  a.  is  the  stress  lever  at  which  dilatation 
is  completely  suppressed;  or,  conversely,  the  stress  level  where  vacuoles 
a-e  first  formed.  Also,  a.  .  is  pressure  dependent, ^defining  the  point  of 
vacuole  formation  over  the  pressure  range,.  Thus,,  a  must  be  a  simple 

function  of  o  That  is,  ^ 

*  • '  cr-  • 

?p  -  f  (1C‘2) 

Equation  (Kr2)  tells;  us  how  to  .evaluate  a  deep  in  the  grain.  The 
critical  .stress  acts  in  the  same  direction  as  the  maximum  principal  stress, 
err  and  it  includes  that  portion,  of  the:  inrier-rbore  pressure  that  is  effective 
deep,  within 'the  grain  and  which  acts  in.  the  same  direction  as  o  .  The 
.fraction,  f,t  of  the  inner-’bore  pressure  that  acts: within  the  grain  can. 'be 
obtained  directly  from  existing sstress;  analyses.  Using  this  fraction  the 
critical  stress  becomes,  on  using  Equation  (14)  of  .the  text, 

o  (F/A)  -  fP.  (K-3) 

cr  cr  i 

where  is  the  inner-bore  pressure 


rJ).  At  this  point, 

requires  large 
a  ,(K-6  to  14) 

sitive  to  relatively  small 


K-3 
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The  appropriate  values  of  ap  for  use  in  Equation  (K-2)  ere  obtained 
from  empirically  obtained!  curve  of  ap  versus  pressure,  the  desired  value 
of  ap,  read  from  the  curve  at  a  pressure  numerically  equal  to  o  ,  the 
latter  being  obtained  from  Equation.  (K-3). 


K-l 

K-2 

K-3 

K-4 

K-5 

K-^6 

K~7 

K-8 

K-9 

K-10 

K-ll 

K-12 

K-13 

K-rU 

K-l  5 


Aerojet  Solid  Propulsion  Company 
Report  1341-26F 

REFERENCES 


Appendix  R 


.'Bills,  K.  W. ,  Jr,.,  Campbell,  D.  M. ,  Sampson,  R,.  C,,  and 
Steele,  R.  D.,  ''Failures  in  Grains  Exposed  to.  Rapid  Changes 
of  Environmental  Temperatures",  Aerojet  Report  1236-81F 
(Contract  No.  N00017-68-C-4415).  (April  1969). 

Bills,  K.  W..,  Jr.,  et.  al,  "Solid  Propellant  Cumulative  Damage 
Program",  Final  Report  No.  AFRPL-TR-63-131 ,  Contract  No. 
F04611-67-C-0102,  (October  1968). 

Wiegand,  J*.  H. ,  "Study  of  Mechanical  Proprccies  of  Solid 
Propellants",  Report  No.  0411-10F,  Aerojet- General  Corporation, 
1S62. 


Hazelton,  I.  G. ,  and  Planck,  R.  W. ,  "Propellant  Characterization 
for  Firing  and  Flight",  Bulletin,  of  45th  Meeting  of.  1CRPG 
Working  Group  on  Mechanical  Behavior.,  p.  287  (Nov.  16-19,  1965). 

Lim,  C.  K. ,  and  Tschoegel,  N.  W. ,  "The  Effect  of  Pressure  on 
the  Mechanical  and  Ultimate  Properties  of  Filled  Elastomers", 
Bulletin  of  the  Eighth  Meeting  of  the  JANNAF  Mechanical  Behavior 
Working  Group,  p;  153  (March  1970), 

Hughes,  D.  S..  et  al,  J.  App,  t/hys .  21^t  294  (1950), 

Masuoka,  S.,  and  Maxwell,  B. ,  J>  Poly.  Sci,  32,,  131  (1958), 

-J.  Res.  Nat.  Bur.  Suds.,  50,  314  (1953).. 

J.  Appl.  Phys.,  30,  337  (1959). 

Bull.  Amer.  Phys.,  Soc.  II,  5,  203  (I960). 

Modern  Plastics,  35,  174  (1957). 

Trans,  Soc-,  Rheol,,  4_,  347  (I960) , 

J.  Chem.  Phys.,  26,  196  (1957). 

Proc.  Roy.  Soc.,  A253,  52  (1959).. 

Farris,  R.  J. ,  "The  Influence  of  Vacuole  Formation  on  the 
Response  and  Failure  of  Filled  Elastomers",  Trans,  of,  the  Soc. 
of  Rheology,  12:2/315-334  (1968). 


APPENDIX  L 

EFFECTS  OF  PREVIOUS  DAMAGE 


Aerojet  Soil’d  Propulsion  Company 
Report  1341-26F 

APPENDIX  ,L 

EFFECTS  OF  PREVIOUS  DAMAGE 


The  cumulative  damage  relations  can  be  used  to  evaluate  the  effects 
of  previous  damage  upon  propellant  failure  data.  This  may  be  done  using 
the  LCD  analysis' of  a  single  test-to-failure  curve  which  duplicates  a  given 
•failure  condition  in.  the  motor.  When  considering  the  inner-bore  failure 
of  a  grain  on  motor  pressurization  we  use  a  high  rate,  uniaxial  tensile  test 
performed  under  a  superimposed  hydrostatic  pressure., 


The  cumulative  damage  relation  appropriate  to  this  analysis  is  derived 
here,  starting  with  Equation  (16)  from  the  text.  For  this  analysis, 
considerations  are  simplified'  by  assuming  the  pressure  and  temperature 
to  be  held  constant.  Thus  a  and  a_  become  constants  arid  may  be  taken 
outside  the  integral  in  Equalion  (16) . 


For  the  previously  undamaged  specimen  the  time-to-break,  t^,  and  the 
maximum  true  stress  become  normalizing  parameters  for  the  Integral  term, 
as  originally  shown  In  Reference  (L-l).,  Thus,  we  let 


A 


* 

o 


o  ) 

cr 


B’ 


d  (t/tfe) 


(Erl) 


Equation  (16)  after  normalization  by.  substituting  A,  gives 


ED  -  1  - 


,^°tM. "  acr' 


CbA 


P(n)  (ato  -  acr)  t^ 


(L-2) 


Prior  to  failure  we  can*  use.  the  same  normalization  terms ,  t^  arid 
(fftM  “  0  r'  ’  'but  take  the  quantity  A  prior  to  failure  with  the  upper 
integration  limit  the  time  (t  <  t^)  at  which  we  wish  to- assess  the 
damage.  Thus, 


t/t. 


A(t)  -  / 
o 


,(c,t  -  ,w 


B 


(atM  ~  °cr\ 


d  (t/tb) 


(E-3) 


Aerojet  Solid  Propulsion  Company 
Report  1341-26F 


Appendix  L 


and 


ED(t), 


°cr)  .  yw 

P<">  <°to  -  acr>B  W, 


(L-4) 


This  relation  is  greatly  simplified  on  dividing  Equation  (L-4)  by 
(L-2)  to  give 


ED(t)  *  A(t)/A 


(I  *5) 


Considering  Equation  (L-5) ,  if  the  specimen  which  had  been  previously 
damaged  fails  at  the  time,  t,  then  the  total  damage  on  the  specimen,  by 
definition,  must  equal  one.  .Since  ED(t)  is  less  than  one,  the  difference 
must  be  accounted  for  by.  the  previous  damage,  ED  .  Thus,  ED,(t)  becomes 


ED(t)  -  1  -  EDp  (L-6) 

The  values  of  A  and  A(t)  can  be  obtained  by  integration  of  the  data 
from  a  single  uniaxial  test  record,  from  which  the  ratio  A(t)/A  and  hence 
ED  can  be  plotted  versus  testing  time,.  From  the  original  test  record  we 
can  make  a  plot  of  the  stress  versus  time  arid  the  strain  versus  the  time. 
Using  these  three  plots  then  for  if  allure  at  any  given  testing  time,t., 
corresponding  values  for  the  stress,  the  strain  and  ED  can  be  read  off » 
Plotting  the  values  of  a  and,  c  versus  ED  gives  the  delired  failure  curves 
and  shows  the  reduction  in  the  undamagedp values  and  caused  by  the 

previous  damage. 
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APPENDIX  M- 

INPUT  DATA  FOR  PRESSURIZATION  TESTS  ON  A  PBAN  PROPELLANT' 

The  data  input  procedures  follow  those  shown  in  Appendix  A.  The 
Prony.  Series  constants  are  tabulated  in  Table  M-l.  The  time-temperature 
and  the  time-pressure  shift  factors  are  tabulated  in  Table  M-2, 

The  cumulative  damage  parameters  for  this  propellant  are  as  fpllows: 


B  =«  8.75 

°to  =  71  psi 

t  *  1  min, 
o 


For  reference  purposes  we  have  provided  a  master  relaxation  cutye, 
Figure  M-l,  and  the  curves  for  a^,  and  ap,  Figures  M-2  and  M-3,  respectively, 


Aerojet  Solid,  Propulsion  Company 
Report  1341t26F' 


Appendix  M 


T’-JLE  M-l 

'  A 

MODULUS-.JLFUT  FOR  PBAN  PROPELLANT 
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TABLE  M-2 

SHIFT  FACTORS  FOR  P BAN  PROPELLANT 


TIME-TEMPERATURE  SHIFT  FACTORS' 


Temp.,.  °F 

lo*10aT 

180 

-2.508 

150 

-1.903 

140 

-0.876 

77 

0 

40 

2.097 

0 

4.057 

-40 

7.008 

-75 

10^528 

TIME-PRESSURE  SHIFT  FACTOR 

Press. , 

gsig-  i 

lo*10?p 

0 

1 

0 

200 

51 

1.71 

600. 

150 

2.176 

1000 

250 

2.398 

-  -  . 

M-3 
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TIKJ5-PRESSURE  SHIFT  FACTOR  FOR  PBAN  PROPELLANT 


— JL _ _ _ . • _ « JL..,.  - _J 
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Figure  M-3 


